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Bell Telephone Laboratories may 
have designed the answer: Nike 
Zeus, a fully automated system de- 
signed to intercept and destroy all 
types of ballistic missiles—not only 
ICBM’s but also IRBM’s launched 
from land, sea or air. The system is 
now under development for the 
Army Ordnance Missile Command. 

Radically new radar techniques 
are being developed for Nike Zeus. 
There will be an acquisition radar 
designed to detect the invading mis- 
sile at great distances. And a dis- 
crimination radar designed to dis- 
tinguish actual warheads from 


harmless decoys that may be in- 
cluded to confuse our defenses. 

The system tracks the ICBM or 
IRBM, then launches and tracks the 
Nike Zeus missile and automatically 
steers it all the way to intercept the 
target. The entire engagement, from 
detection to destruction, would take 
place within minutes and would span 
hundreds of miles. 

Under a prime Army Ordnance 
contract with’ the Western Electric 
Company, Bell Laboratories is 
charged with the development of the 
entire Nike Zeus system, with assist- 
ance from many subcontractors. It 
is another example of the coopera- 
tion between Bell Laboratories and 
Western Electric for the defense of 
America. 
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STATISTICAL DESIGN OF DISCRETE-DATA CONTROL 
SYSTEMS SUBJECT TO POWER LIMITATION 


BY 


JULIUS T. TOU' AND K. S. PRASANNA KUMAR! 


SYNOPSIS 


A solution is presented for the analytical design of a digital controller in ac- 
cordance with a specified optimum criterion. This approach makes use of the 
z-transform and the modified z-transform techniques. The analytic design procedure 
for discrete-data control systems subject to power limitation is developed. Attention 


is centered upon stationary random functions because of their frequent occurrence 


in control systems and their being amenable to mathematical treatment. In this 


paper, physical realizability and system constraints are considered and general 


design equations are derived. 


INTRODUCTION 


Conventional design techniques where the system inputs are as- 
sumed to be deterministic in nature may fail to yield satisfactory re- 
sults. This is due to the fact that the designer presupposes the nature 
of the input as either sinusoids or simple aperiodic functions. In order 
to broaden the design theory, stochastic input signals should be con- 
sidered (1, 8, 9).2 

The analysis and synthesis of discrete-data control systems has 
been a matter of considerable interest for well over ten years (1-7). 
Several papers have appeared on the statistical design of digital and 
sampled-data control systems. These papers (1-3) have formulated 
systematic procedures for the design of discrete-data systems satisfying 
a minimum mean-square-error criterion. This paper extends the con- 
cepts and ideas of the author’s previous work (1-3) to the design of a 
discrete-data control system such that the mean-square-error is mini- 


1 The Information Technology Laboratory, the Technological Institute, Northwestern 


University, Evanston, III. 
2 The boldface numbers in parentheses refer to the references appended to this paper. 


(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
the JouRNAL.) 
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mized while at the same time the power of the control signal is con- 
strained. This paper may be considered as the companion paper to 
reference 3 which appeared in this JOURNAL. 

Power constraint is a very important aspect in control-system design. 
Without considering power limitations, the design may lead to a system 
with excessively large power dissipation in the components (8, 9), 
which is certainly undesirable and impractical. This paper is con- 
cerned with this important design factor. Optimum pulse-transfer 
functions for the digital controller are derived under physical realizabil- 
ity conditions. The input signals and noise are assumed to be sta- 
tionary random functions. 


OPTIMUM DIGITAL CONTROLLER FOR MINIMUM MEAN-SQUARE SAMPLED ERROR 


The optimum design problem is solved with reference to Fig. 1. 
The mean-square sampled error of the system is to be minimized, 
meanwhile the mean-square value of the control signal u(”7°) is con- 
strained. Now, let G.(z) be the over-all pulse-transfer function of the 


e(nT) 


Fic. 1. 


system which is to be optimized in the above sense. (Ga(z) is the de- 
sired over-all pulse-transfer function. The sampled input and output 
of the system are: 


r(nT) = r,(nT) + 7,(nT) (1) 


and 


c(nT) = c,(nT) + ¢,(nT), (2) 


respectively. 
The desired output sequence is ca(n7). The error sequence of the 
system is given by 


(3) 


e(n7T) = c(nT) — ca(nT). 


The mean-square sampled error is given by 


——— 1 
e(nT) = af, dee (4) 


+ 
Weg 
fal 
r, (nT) C,(nT) 
: 


Sept., 
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in which 
= [Go(2) — Gal) — Gale) 


+ G.(2")[Go(z) — Ga(Z) + (5) 


value of the signal u(”7) is given by 


1 
u2(nT) burn (6) 
in which 


duu(3) = (3) Drarn(S) + Prern(Z) + Prarg(Z) (7) 


The mean-square 


Referring to Fig. 1, 


= W(z)G(z) (8) 
where 
D(z) 
Ws 9 


and, G(z) has its poles and zeros inside the unit circle. Substituting 


(8) into (5) yields 
dee(s) = [W(2)G(s) — Gals) — Gals) (2) 

+ — Ga(s) (3) 

+ W(s7)G(s") LW (2)G(z) — Gals) (2) 

+ (2). (10) 


The design objective is to minimize (4) subject to the constraint 
that (6) is less than or equal to a specified value. u?(77) may repre- 
sent the power input to a component of the system and it may be 
necessary to constrain it for protection of the part concerned. 


The minimization problem is solved as follows. The functional 


F=e(nT) + hu? (nT) (11) 


is formed where 


\ = the Lagrange multiplier. (12) 


Substituting (4) and (6) into (11) yields 


1 
Jr 


2 

i 

‘ 
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Following the method of reference 1, the pulse-transfer function W(z) 
is given a small variation yn(z), and W(z~') is given a variation yn(z~). 
As a result, F undergoes a variation 6F. (zg) is an arbitrary function 
of gand y isa parameter. W(z) is replaced by its neighboring function 
W(z) + yn(z), and W(z-') by its neighboring function W(z-) + yn(z-'). 
The functional F becomes F + 6F. Then, F is a minimum, if 


_ 
(14) 


Following the operations involved in (14) and simplifying leads to 


n(z){W(z—)[A + G(s)G(2-) ]6(z) — G(s)Ga(s—) (2) 
Jr 


+ (2) {W(2)[A + 
— = (15) 


where 
(2) A + Prern(Z) + brare(Z) + brarn (3) (16) 
$1(2) A + (3) (17) 
A + or, r, (18) 


The next step is to factor ¢(z) and [A + G(z)G(z~!) ] into two parts: 
$(z) A (19) 
[A + G(z)G(z-') ] 4 B(s) = Bt(z)B-(s), (20) 


where the first factors of (19) and (20) contain zeros and poles inside 
the unit circle and the second factors contain zeros and poles outside 
the unit circle. From (15), (19) and (20) results 


(2) Bt(z) 


= 0. (21) 


(z)B-(z) 


ij 
‘ 
3 
: 
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For physical realizability, W(z) and n(z) can have no poles outside the 
unit circle and W(z-') and n(z~!) can have no poles inside the unit circle. 
Thus, the function 7(z)¢+(z)B+(z) contains poles inside the unit circle 
only, and the function n(z~')@~(z)B-(z) contains poles outside the unit 
circle only. Consequently, (21) may be written as 


n(z)p* (2) Bt (z) (2) B-(z) 


ot(z)Bt(z) J- 


mJr 


| B-(2) idg=(0 (22) 


in which the symbol { }, means the operation of picking the part of a 
function of z with poles inside the unit circle, and the symbol {| }_ means 
the operation of picking the part of a function of z with poles outside 
the unit circle. For instance, in F(z) = {F(z)}, + {F(s)}_, the terms 
with poles inside the unit circle are grouped in the (+) part, and the 
terms with poles outside the unit circle are collected in the (—) part. 

If W(s) is indeed the optimum transfer function, (22) must be 


satisfied for arbitrary n(z). Thus, the optimum W(z) is given by 


$-(z)B-(z) 
(2) Bt (z) 


(23) 


W(z) = 


The pulse-transfer function D(z) is determined from 


W (sz) 


(24) 


D(z) = i 


It is noted that (23) and (24) contain \ asa parameter. The value 
of \ is adjusted so as to satisfy the constraint on u?(n7), namely, 


u2(nT) < (25) 


In practice, the adjustment of \ may entail a graphical procedure if the 
constraint equation cannot be solved explicitly for \. When once the 
value of \ is fixed, (23) and (24) yield the optimum pulse-transfer 
functions. From these and (4), the mean-square sample error can be 
calculated. 


: 
She 
| 
: 


176 


Juttus T. Tou anp K. S. Prasanna KuMAR (J. F. 1. 


Example 


1 Ts 
Given: G(s) = re Ga(s) = 1; T = 0.2 sec. 


rer, (@) 


(w) 


u2(nT) < 2.0. 


Determine the transfer function of the pulsed-data compensator re- 
quired to minimize the mean-square sampled error. The following 
s-transforms follow. 


0.181 


—(), 20132 
— 0.819) — 1. 221)’ 


— 0.261) (2— 3. 


Prer, (2) 


— 


~ 


From (20), 


(2 — a)(s — b) 
s — 0. 819) - 1.221) 


{ (2.04 + 2%) {(2.04 + 2%)? ~4l 
{ (2.08 + 2%) + (2.08 + 


when \ > 0, a < 1 and b > 1 and hence 


B(z) = 


where 


B+(z) = 819) and B-(z) = 1.221)" 


Using (23), 


0.1(s — 0.819) 


— 0.819) (z — a)(z — 0.261)" 


: 
1 
Prara(w) = 0.1 
3 
2412 
24 
= 0.1 
791) 
791 
221) 
: 
and 
b= 
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From (6), 


—— 0.00314 (3.791 — 
w(nl) = 


anr?(b — 0.819)? (- (3.83 — a) 


Actually a curve of u2(nT) v/s e?(nT) with as the parameter should 
be plotted and the optimum X selected therefrom. However, if \ = 0.01, 
u?(n7T) = 1.25 which satisfies the specification u?(n7T) < 2.0. Thus, 


1.98(21 — 
(1 — 0.4372-' + 0.0442-*)’ 


1.98 (s-' — 0.819z-*) 
(1 — 0.4372—! — 0.314s-2)’ 


“a (nT) e(t)=C (t)+c, (t) 
r(t)=r,(t) Sta) n n 


G,(s) 


Fic. 2. 


OPTIMUM DIGITAL CONTROLLER FOR MINIMUM MEAN-SQUARE ERROR 


With reference to Fig. 2, 


e(t) 


|! 


c(t) — ca(t) 
c,(t) +- c,(t) ca(t). 


(26) 


In terms of sampled values, 


(27) 


e(nT, m) = c,(nT, m) + ¢,(nT, m) — ca(nT, m) 


where 7 is an integer and m lies between 0 and 1. For a given value of 


m, the mean-square value of the error samples is 


: 
177 ‘ 
fe 
, 
= 
J ( 2) 
from (9) that 
It follows from ( 
: 
D(z) = 
- + e(t) 
: 
C.(t 
t 
r ( ) 
s 
= 
4 
= 
: 
: 
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e (nT, m) = ca(nT, m) ¢,2(nT, m) + c,?(nT, m) 


+ m) + ca(nT, myc.(nT, m) 


— ca(nT, m)c,(nT, m) — c.(nT, m)ca(nT, m) 


— ca(nT, m)c,(nT, m) — c,(nT, m)ca (nT, m). (28) 


Equation 28 reduces to 


e(nT,m) = dee (Zz, m)z—dz (29) 
T 


where 


Pee(Z, M) = Peaca(Z,M) + Pegcg(Z, M) + Henen(Z, M) 
+ (2; m) Penes (z, m) =4 Peacs (z, m) Pesca m) 
— eac,(Z, M) — M). 


(30) 


If G.(z, m) is the over-all modified pulse-transfer function of the system 


G.(s) = W*(s)G(s) (31) 


where 


W*(s) = D*(s)/[1 + D*(s)G*(s) ] (32) 


we have 


G.(z, m) W(z)G(z, m). (33) 


The block diagram of Fig. 2 is redrawn in Fig. 3 with the feedback 


L 
r(t)er (t)+r, (t) | 
T 


q 


Fic. 3. 


system represented by its open-loop equivalent and the desired transfer 
function G(s) preceded by a sampler and an ideal desampler. In order 
for the signal input to Ga(s) to have the same statistical properties, it 


: 
| 
G (s 
(s) 
u(nT) 
W (s) t G(s) C(t) 
+ 
e(t) 
O 
- 
r(t) (t) 
s H s t 
s (s) G 4 (s) 
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can be shown that the transfer function //,(s) of the ideal desampler 
is given by 


H, (s) (s). (34) 


In view of this, (30) reduces to 


dee(z, m) = m)G(z-, m) 
— W(z)G(z, m)Ga’ m) 
— W(z-)G(z-!, (z-', m) 
— (2, m)Ga' (2, Mm) (2) 
+ [W(z)W(z—)G(z, m)G(z—, m) 
— W(z)G(z, m)Ga' (2+, m) 
+ [W(s)W(s—)G(z, m)G(z-1, m) 
— W(s")G(s-, m)Ga' (2, m) (2) 
+ m)G(z-, m)¢r,+, (2) 


(35) 


where G,'(z,m) is the modified z-transform associated with G,’(s) 
= H,(s)Ga(s) and (Z), Prer,(Z) and ¢,,,r,(Z) denote the pulse- 
spectral densities and the cross pulse-spectral densities for the input 
signal 7,(¢) and the input noise 7, (t). 

From Fig. 2, 


1 
u(nT) = (36) 


where 


duu(s) = + Prarn(Z) + + (37) 


Now form the functional 


F = e&(t) + w(nT) (38) 
where 
e?(t) -f e?(n7T, m)dm (39) 
1 td f 1 
J dee(Z, mdm. (40) 
We have 


— — 
+ (2) 


: 
Gps 
i 
: 
(4 1 ) 
i 
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where 


G(z, m)G(z-', m)dm 


K,(z) 


K,(z) (zs, m)Ga'(s~', m)dm 
K;(z) f 
4+ 


Ga’ (z, m)Ga’ m)dm 


= 


+ + (8) 


= (2) + (2) 


and 


The minimization of (38) subject to the constraint that 
w(nT) < 0? 
is solved as before. The result is 


K2(z Lor, ra ) Ts, (z) J | 
(z 
W (sz) 
where B(z) = [A + Ki(s) ]. 


SYSTEMS WITH NOISE IN THE CONTROL LOOP 


With reference to Fig. 4, the system error is given by 


ca(t) — c(t) 
Ca(t) — c,(t) — cn (t). 


(43) 
(44) 
3 (45) 
6 
(46) 
(47) 
: 
(48) 
(49) 
« 
= (50) 
n(t) 
D*(s) 
r,(t) u(t) “u(t) e(t) 
: 
e(t) 
O 
: r,(t) C(t) 
: 
; Fic. 4. 
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In terms of the error samples, 


e(nT,m) = ca(nT, m) — c,(nT, m) — c,(nT, m). (51) 


The case where the input signal 7,(¢) and the noise n(¢) are uncorrelated 
will be considered. The mean-square value of the error samples is 


e(nT, m) = ca(nT, m) +.,2(nT, m) + ¢,2(nT, m) 


— c,(nT, m)ca(nT, m) — ca(nT, m)c,(nT, m). 


(52) 


The mean-square value of the error samples is given by (29), where 


= hegeg(Z, M) + + Hence, (Z, Mm) 
— M) — deac,(Z, Mm). 


m) 


(53) 


Figure 4 is redrawn in Fig. 5 to simplify the evaluation of ¢..(z, m), 
where /7,(s) and /7,(s) are the transfer functions of the ideal desamplers 


(s) 


W"(s) G(s) 


G(s) 


Ho(s) 


Fic. 5. 


The transfer function H,(s) 


for the signal 7,*(t) and the noise n*(t). 
is given by (34) and it can be shown that 


H,,(s) = (54) 


8 
Io! 
4 
: 
: 
: 
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The modified pulse-spectral densities of (53) are given by 


(2,m) = W(s)W(s-)G(z, m)G(z-', m)¢,,r, (2) (55) 


Pegcs 


deaca(Z, M) 


= Ga’ (z, m)Ga' (2) (56) 


Peyca(Z, m) = m)G(s-', m)Ga' (z, m)d,,r, (2) 


Deac,(Z,m) = W(s, m)G(z, m)Ga’ m)d,,r, (2) (58) 
W (z)W (2-)G(z, m)G(2-, m) dain, (2) 
— W(z)G(z, m)G2' m) dan, (2) 
— W(s)G(z-', m)G2' m)dnin(2) 
+ (s, mM)dan(z). (59) 


PDenen (z, m) 


In (59), G.'(z, m) is the modified z-transform associated with 


(s) = = (s)\Go( s) (60) 


and 


Ni(s) = €n,(s). (61) 


In Fig. 4, 


G(s) Gi(s)Go(s). (62) 


Referring to Fig. 4, the constraint is 


u2(t) < (63a) 


which can be written as 


u(t) < (630) 
where 


~ 

= | 
t 


= o” 2(t) (63c) 


and 


u;(t) = = 


/ 
du (2, m)z—'dz. (64) 


In (64), du,u,(Z, m) is given by 


Pusu m) = m)G,(s-, m)¢,,,, (2). (65) 


Application of the calculus of variations to the functional F = e*(t) 
+ u,°(t) yields, 


B-(z 
W(s) = (66) 


(57) 
i 
Berd 
is 
> 
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where 


B(s) = + Gnn(Z)] + (2) (67) 


K,(z) = G(z, m)G2'(z—, m)dm (68) 


and 


K,(z) = f G,(z, m)G,(z-, m)dm. (69) 


Clearly, the technique could be extended to cases where constraints are 
placed upon more than one signal in the control loop and the signal and 
noise are correlated. 


CONCLUSION 


This paper introduces a simplified systematic procedure for sta- 
tistical design of discrete-data control systems when the control signal 
is constrained. The system performance is optimized in the sense of 
minimum mean-square error. The system inputs are assumed to be 
functions of stationary random processes. Although only systems with 
one constrained signal are discussed, the design technique could be 
easily extended to cases where constraints are placed upon more than 
one signal in the control loop. The application of this procedure to 
optimum design for minimum integral-square error and subject to 
power limitation is apparent, if use is made of the results obtained in 
reference 3. Furthermore, this design technique may be extended to 
certain types of nonlinear control systems, provided that the non- 
linearity can be described or approximated as constraints on control 


signals. 
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STRESSES IN CYLINDRICALLY SYMMETRIC MEMBRANES 
REINFORCED WITH EXTENSIBLE CORDS 


BY 


W. F. AMES! 


SUMMARY 


The theories of Rivlin and,Hofferberth on the equilibrium shape and stresses in 
deformed membranes reinforced by cords are limited to inextensible cords and to an 
average cord-extension ratio, respectively. The results of this paper remove these 
limitations for the case of cylindrically symmetrical deformation. The equilibrium 
shape and the stress problems are solved for the deformed membrane. 


1. INTRODUCTION 


In recent years there have been a number of quite successful efforts 
to apply the membrane and shell theories of elasticity to the analysis 
of the deformation and stresses in reinforced structures. Of special 
interest in this area is the work of Adkins (1, 2)? Rivlin (3, 4), Adkins 
and Rivlin (5), Hofferberth (6, 7, 8) and Lauterbach and Ames (9). 
In their earlier papers (1, 2, 5), Adkins and Rivlin considered the de- 
formation of sheets and shells of a material consisting of incompressible 
isotropic elastic material reinforced by inextensible cords. The pres- 
ence of the elastic material made it difficult to solve any but the simplest 
problems. 

A number of technological problems, of interest to engineers, are 
such that either the elastic material is absent altogether (for example, 
cloth, canvas, etc.) or the elastic material plays a subordinate role in 
the determination of the deformation and stress characteristics of the 
system. Examples of the latter case include the pneumatic tire, bal- 
loons, radomes and so forth. Realization of this subordinate position 
of the elastic material prompted Rivlin (3) to consider the plain strain 
of a net formed by two families of straight parallel inextensible cords. 
As an outgrowth of this work a fundamental paper by Rivlin (4) con- 
sidered the deformation of a membrane formed by inextensible cords. 
A direct application to engineering structures is possible from this 
theory, especially that case of cylindrically symmetric deformation 
which concerns us here. 

Paralleling, but independent of, the work of Rivlin, were the efforts 
of Hofferberth (6, 7, 8). Hofferberth was entirely concerned with the 
development of an analytic tire theory based upon membrane theory. 
Consequently, his results are not as general as those of Rivlin. Early 


1 Department of Mechanical Engineering, University of Delaware, Newark, Del. 
2 The boldface numbers in parentheses refer to the references appended to this paper. 
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in this work Hofferberth’s results, on the membrane equilibrium shape, 
became known to the present author. By use of these results the stress 
problem in cylindrically symmetric deformation (tirelike structures) 
was solved independent of the work of Rivlin. Lauterbach and Ames 
(9) demonstrated that the tire stresses in the two principal directions 
and the cord tension are quite variable which results from the expan- 
sion process. 

Also to be mentioned in this area of work is the book of the Russian 
V. L. Biederman on ‘Tire Design.” 


UNDEFORMED SHAPE 
DEFORMED SHAPE 


Fic. 1. 


The Rivlin theory has the defect that the cords are assumed in- 
extensible. Hofferberth’s theory purports to remove this limitation but 
his membrane shape includes a cord-extension ratio whose dependence 
on the horizontal dimension y is unknown. Thus his results are limited 
to the case of using an average cord-extension ratio. The departure 
from reality, in this case, is not severe since the cord strain realized in 
tires is of the order of 2-4 per cent resulting from standard inflation 
pressures. However, in other applications, such as reinforced tubing, 
this limitation may be extremely important. 

The present paper is concerned with removing the inextensible 
limitation of previous work for the case of cylindrically symmetric 
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deformation. A Hookean relation is assumed between the cord strain 
and the cord stress. The cord extensibility is then obtained in de- 
pendence upon dimension and the equation for the shape of the de- 
formed membrane is derived. Finally the stress problem is solved 
illustrating its dependence upon physical parameters, membrane shape 
and type of loading. 


2. BASIC EQUATIONS 


The notation used throughout the paper is consistent with that of 
the previous paper (9) which agrees with that of Hofferberth (8). 

The deformed membrane is a surface of revolution so that an or- 
thogonal net with coordinates (¢, 6) is chosen on the surface. Rec- 
tangular coordinates (y,z) are chosen in the plane as illustrated in 
Fig. 1. We shall refer to @ as a meridian coordinate and @ as a circum- 
ferential coordinate. 

From (4, 8, 9) the meridian equation is 


dz 
dy = (1) 


and the shape equilibrium equation of the membrane is 


do 4 cot? B iv 2y (2) 


‘ot 


where 8 = B(y) is the angle of the deformed membrane cords with the 
circumferential circles and a is the constant cord angle with the circum- 
ference of the membrane in the predeformed cylindrical shape. 

From the cord network considerations of Rivlin (3) as modified by 
Hofferberth (8) we have the relation 


Ci, COS B (3) 


y= 
COs @ 


where c, is the cord-extension ratio and the other parameters are defined 
in Fig. 1. 
From (9) the cord tension r is expressible in the form 


(y? — cot (4) 
sin sin a cos a 


where p = internal pressure, 27) = number of cords per inch of material 
before deformation and N = number of cord layers. 
Finally, from the membrane theory of elasticity, the equations of 
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force equilibrium for the membrane are 


10(Ty) ON, 

+ = + T cos + pey 
10(Noy) , OT 
al 0d 

a) 


No cos -}- Poy = 


where N, = stress component in the circumferential direction, 
N, = stress component in the meridian direction, T = Nye = Neg 
= shear component, pp = load component in the circumferential direc- 
tion, p, = load component in the meridian direction, p, = load com- 
ponent normal to the surface, and 7, 72 = radii of curvature in the 
meridian and parallel direction, respectively. 


3. CORD EXTENSION DEPENDENCY ON y 


We now assume a linear relation between the cord stress 7 and cord 
strain e, € = Az so that the cord extension ratio can be written as 


¢.=1+e=1+Ar (6) 


where A is the appropriate physical constant. With the assumption 
of (6) the occurrence of c, in the shape integral (see 8, 9) can now be 
eliminated by obtaining the dependence of ¢; on y. 

From (3) 


y COS @ 


=1+Ar. (7) 


Cy = 
r, cos B 


By using (4) and (7) there results, after some algebra, 


1 E COS a (y? — 


AlLr, cos B 2Nnor, sin B sin ¢ sina 


Solving for sin ¢, we have the result 


(y? — COS 
sin B(y cos a@ — cos B) 


where 
p 
(10 
2Nn of, sin aA 
and (9) is independent of c, but does still depend upon 8. The de- 
pendence of 8 upon y must be found in order that sin @ and c, will have 
an exhibited relationship with y. 
Upon differentiation of (9) with respect to y and division by sin 
there results 


. 
| (J. F. 1. 
= 
(S) 
“3 
: = 
(8) 
: 
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cosa+r 
cot — (tan B + cot B) ay, (11) 


dy y— yz ly ycosa—~r,cosB 


But, substituting (2) in (11) there results 


d 
cosa + 7, sin B- 
cot? B dy 
— = (tan B + cot “> — (12) 
y ycosa — r,cos B 


which is the differential equation for B as a function of y. The solution 
of (12) is facilitated by the introduction of the new variables 


u =cosBp (13) 


y cos @ 
rs 


whereupon (12) becomes 


du ul v(2u? — 1) — u*] 
— v) 


which has the solution 


= 1). 
(15) 


When y = ¢,, no expansion of the cord matrix has taken place from the 
undeformed cylindrical state. Therefore, at y = 7,, 8 = a and u, = 2, 
= cosa. It is easily seen that C = 3 and (15) becomes 


— (2v? + 1)u? ++ = 0 (16) 


which is the implicit expression for u = cos B as a function of v, and 
hence of y. Clearly 0 < u <1 so that the desired root of (16) lies 
between 0 and 1. The function given by the left-hand side of (16) is 
easily shown to have a minimum (of interest here) at u = —— with 


the minimum value 


* 


Returning to (9) we write 


(1 — — u) (17) 


sin @ = 


pr. Vp COS 
and vg = ——. As v varies from 


2A Nny sin a cos? a@ Ys 


where B = 


2 
18 
9 
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cos a to 7, cos a/r,, the corresponding value of « = cos 8B may be found 
by numerical solution of the cubic equation (16). Hence c, may be 
found from (7). 

When (17) is substituted into (1) the shape of the deformed mem- 
brane becomes 


r,B(v? — vpz*)udv 
cos al (1 — u?)(v — u)? — B2(v? — 


where the negative sign applies when z is positive and the positive sign 
when gz is negative. The quantity 1. = — is introduced as the 
& 
lower limit since z(7,) = 0. This integral is hyperelliptic and evalua- 
tion numerically is necessary. 


4. MEMBRANE STRESSES 


The solution of the membrane equilibrium equations (5) has al- 
ready been given in general form in (9). If we denote by f(y) the func- 
tion under the integral sign in (18), that is f(y) = tan ¢, we have from 
(9) that the stress N, in the circumferential direction is 

df 


(y? Vs") 


dy (19) 


_ pyV1 + f(y) 
f(y) 2y f(y) (1 + f2(y)) 


Neo 
and the stress N, in the meridian direction is given by 


FPO) 
f(y) 


Finally, the tension 7 in the cords as a function of y is given by 


p(y? — + f?(y) 


Nnor, sinasin8  2Nnor, sin af(y)¥1 — cos? 
where cos 8 is calculated from (16). 
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ANALYTICAL SOLUTIONS FOR ROTARY MATRIX 
WIRE SCREEN HEAT EXCHANGERS 


BY 
W. H. T. LOH! 


ABSTRACT 


Instead of the partial differential equation commonly used for rotary matrix heat 
exchangers, an ordinary differential equation is developed for a single stage rotary 
wire screen heat exchanger. This equation results in an exact solution for tem- 
perature effectiveness. The solution for a multi-screen system is then developed 
by an energy balance of each individual single screen. Based on analytical solutions 
obtained here, design charts for rotary multilayer wire screen heat exchangers are 


presented. 


INTRODUCTION 


The rotary matrix heat exchanger, developed in early 1922, em- 
bodies a moving heat conveyor, heated in the hot stream and cooled in 
the cold stream. Heat is thus transferred by rotating the conveyor 
through the two streams. The heat conveyor originally proposed was 
a rotary matrix made of plain metal plate surfaces. Improvements in 
heat transfer were effected through the use of corrugated plate-fin type 
surfaces. Further improvements may be obtained by using multi- 
layer individual wire screen as the conveyor. 

Analytical studies relating to rotating heat exchangers have been 
made by Nussett, Hausen, Boestad, Saunders and many others. 


HOT GAS FLOW 


COLD GAS FLOW 


Fic. 1. Flow passage length for generalized rotary matrix. 


However, these authors have concentrated on the solution of the partial 
differential equation which was derived primarily for the generalized 
rotary matrix along the lines of the original proposed arrangements. 
This type of arrangement has a flow passage length (ZL) as shown in 
Fig. 1. 

It can be seen that the temperature (7) is a function of both time (t) 
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and space (x), and that the only equation describing this type of matrix 
is partial differential equation 


oT aT OT, hA 


= — f,). (1) 


Because of mathematical complexities, solutions were generally ob- 
tained either by graphical or numerical methods, since no analytical 
closed form solution has been found. The present analysis takes a 
different approach. This approach, primarily for multilayer rotary 
matrix, starts with a single screen (see Fig. 2). 


HOT GAS FLOW 


COLD GAS FLOW 


— SINGLE SCREEN 


Fic. 2. Flow passage length for single screen 


From Fig. 2, it can be seen that the temperature (7°) is a function of 
time (¢) only. Therefore, an ordinary differential equation was de- 
veloped, resulting in an exact and closed form solution. The multi- 
layer screen system was then developed by an energy balance of each 
individual single screen. Since this approach led to a closed form 
analytical solution for the rotating multilayer wire screen heat ex- 
changer, it was then possible to develop design charts for this type of 
heat exchangers. 


ASSUMPTIONS 


1. Gas and material properties such as p, C>, k, uw, P,, etc. are treated 
as constants. 

2. Heat transfer through conduction and radiation is small com- 
pared with heat transfer through convection ; conduction and radiation 
heat transfer is, therefore, neglected. 

3. Variations in the gas temperature ahead of each individual screen 
are small when compared with the gas temperature itself, hence the gas 
temperature ahead of each screen is treated as uniform, as is the gas 
temperature after each screen. 

4. The size of wire is so small that it may be treated as a point for 
gas flow across it. 
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ANALYSIS 

Equal Mass Flow for Hot and Cold Streams 
Equal mass flow in equal cross-sectional areas results in equal values 

of G, and consequently, equal values of h in both hot and cold streams. 
1. Wire Temperature vs. Time 


For a cross-wire screen element as shown in Fig. 3, the energy equa- 


tion gives 


dT, 4h ) — 4h )r 
(2) 


Fic. 3. Cross-wire screen element. 


As the wire element rotates from hot stream to cold stream, the sur- 
rounding gas temperature (7), shown in Fig. 4 can be expressed by the 
Fourier series 


T= - sin nwt + . (3) 


n=odd 


Solving Eqs. 2 and 3 gives 


Fe = bo + Qn F + (nw/c)? 


n=odd 


n=odd V1 + (nw/c)? 


Equation 4 describes the complete time history of the wire temperature 
from the very start to steady state cyclic operations. Since the tran- 
sient portion becomes negligibly small after an elapsed time approxi- 
mately equal to (3/C) seconds from the very start of the operation, 
only the steady state portion of the complete solution is really needed. 
This steady state solution in non-dimensional form is 


Tw — bo _ > (nwt — dn). (5) 


n=odd nvi + (nw/c)? 
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Equation 5 describes the complete wire element steady state cyclic 
temperature variations in terms of wire geometry, wire material proper- 
ties, hot and cold stream temperatures, heat transfer coefficients, and 
revolutions per second of the wire screen rotor. 


FOR EQUAL MASS FLOW, 


t 


Fic. 4. Gas temperature in rotating wire element. 


2. Cyclic Wire Temperature Rise 


After steady state heat transfer is attained, the zero time (¢ = 0) 
can be set at the instant that the wire element enters the hot stream. 
The wire element leaves the hot stream at t = t,. If 7.,, indicates the 
wire element temperature in °F. at ¢ = 0 and T,, indicates the wire 
element temperature in °F. at ¢ = ¢,, then from Eq. 5 the values of 
T., and T,,, can be calculated as: 


= Mn sin (—¢n) Kn sin Pn 


n=odd n=odd 


sin (nt — = Mn SIN by 


n=odd n=odd 


= tan“'{ — }. 
nv1 4 + /c)? c 


Consequently, the wire temperature rise A7,, from 7,, to 7, during 
the passage of the wire element through the hot stream section is 


2un SiN dy. (8) 


ae 

T, T, 

h Th 

T 

: 

h Cc 

a 

— (6) 

7. 

2 

Mn = 
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3. Temperature Effectiveness of a Single Rotating Wire Screen 


The temperature effectiveness 7, of a single rotating wire screen is 
defined as 


(9) 


Using the energy equation, one obtains 


— Two) = — 


From geometry and the equation of continuity, one obtains 


Substituting Eqs. 8, 10, 11, 12, and 13 into Eq. 9 and simplifying, one 


W, 


obtains 
Here 


¢ may be called the “heat transfer geometry parameter,’’ and 8B may 
be called the “gas flow geometry parameter.”” These two nondimen- 
sional parameters alone determine the temperature effectiveness of the 
single layer wire screen. It should be noted that, when the gas flow 
geometry parameter, 8, is small or approaches zero, the temperature 
effectiveness approaches 1.184¢. 


4. Temperature Effectiveness of a Multilayer Rotating Wire Screen 
Heat Exchanger 


For a single layer wire screen rotor, the results of the expression of 
temperature effectiveness are the same for both parallel flow and counter 


e 
: 
2 
: 
(10) 
: 
¢ 
: 
: 
A 
ae 
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flow heat exchangers. However, the results for parallel flow and 
counter flow are different for multilayer rotating wire screens. These 
results are derived separately as follows: 


a. Parallel Flow 


Figure 5 illustrates a typical arrangement of a multilayer rotating 
wire screen heat exchanger having 7 individual screens. 


HOT STREAM FLOW —> Th 


COLO STREAM FLOW —> 


SCREEN NO.3 
SCREEN NO.4 
SCREEN NO.5 
SCREEN NO.n 


z 
” ” 


Fic. 5. Typical multilayer rotating wire screen heat exchanger having 
n individual screens, parallel flow. 


By definition, 


therefore, 


1 = 


T, 


3 

h 

| hal | : 

n+l 

7 h 1 7 ci 7 h 1 — 7 

Ti, — Th; T,, Th; 

T,, — Ti, T,,- 
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but, 


Te, = Te, + (Ta, — Te). (21) 


Substituting Eqs. 20 and 21 into Eq. 19, one obtains 


n2 =m + — 2m). (22) 

Similarly, 
ns = + — (23) 
Mn = + — 200-1). (24) 


te. = Me, Eqs. 18, 22, 23, 


For the special case of 7., = m, = ™; = 
and 24 reduce to 


Ns 
m + — 2m) 
= + — 2n2) 


Nn-1 + ns(1 2nn—1) 


HOT 
STREAM 


n-2 


L 


SCREEN NO. 4 
SCREEN NO.nN 


SCREEN NO.3 


SCREEN NO.2 


Fic. 6. Typical multilayer rotating wire screen heat exchanger having 
n individual screens, counter flow. 


b. Counter Flow 


Figure 6 illustrates a typical arrangement of a multilayer rotating 
wire screen heat exchanger having m screens. By definition, 


; 
: 
: 
t 
3 
m1 
Ne 
- 
(2 ) 
: 
13 
& | 
: 
COLD 
FLOW 
‘ 
z 
z 
Ww 
: 
: 


Therefore, 


= T,, — T,,) 


Using Eqs. 29 and 30, it is easy to prove that 


Substituting Eqs. 31 and 32 into Eq. 28 results in 


1 — 1 
1 — nem 1 — 


Similarly, 


(34) 


— Ns3N2 Ns 372 


Nn = Nn- — ) + m, (35) 
990-1 1 — 


For the special case of 9., = n.. = 1; = 1, = Me, Eqs. 27, 33, 34 and 
35 reduce to 
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(26) 
= (27) 
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= 1 ( (28) 
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(30) 
(31) 
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ni = Ne 
1 — 
nan 5") 
1 — mm 1 — nm 
1— 8 1— 2 
(2%) ( 
1 — 1 — 
1 —», 1— 
on (753) 
1 — 1 — 


UNEQUAL MASS FLOW FOR HOT AND COLD STREAMS 


| 


3 
| 


(36) 


3 
| 


For unequal mass flows, the same value of G in both hot and cold 
streams may be obtained by designing the ratio of the cross-sectional 
area of one stream to that of the other so that it is equal to the corre- 
sponding mass flow ratio of the two streams. This results in the same 
values of h in both hot and cold streams. 


1. Wire Temperature vs. Time 


Equation 2 is also valid for unequal mass flows. The surrounding 
gas temperature (7), shown in Fig. 7, however, becomes 


= + Ya,’ sinnot + cos nut 


n=l 


T, — 


(1 — cos nbr) 


The solution of Eqs. 2 and 37 is 


Te = + bo + |* (nut + Yn — dn). (38) 
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— 
V1 — cos nbr 


sin nl 
Pn = tan7! ( Wn = sin now | (39) 


1 — cos nbr 


The value of the arbitrary constant A,’ can be determined by setting 
t = Oand 7, = T,, in Eq. 38: 


V1 + (w/c)? 


n=1 


FOR UNEQUAL MASS FLOW, 


t 


Fic. 7. 


After an elapsed time equal to (3/C) seconds from the very start of 
the operation, the transient state portion becomes negligibly small. 
Therefore, for practical purposes, only the steady state portion is 
needed. This solution in nondimensional form is 


: || sin (nwt+y,—¢,). (41) 
| tL nvV1+(nw V2 2 | 


Equation (41) reduces to Eq. (5) when 6 = 1, for equal mass flows. 


2. Cyclic Wire Temperature Rise 


By employing a procedure similar to the one described for the 
corresponding portion of equal mass flows, one obtains the wire tem- 
perature rise: 


‘ 
Here 
: 
7, T 
h h 2 
: 
t 
T 
| T Te 
t 
A 
3 
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kor equal mass flows, Eq. 42 reduces to Eq. 8 when 6 = 


3. Temperature Effectiveness of a Single Layer Rotating Wire Screen 


By employing a procedure similar to the one described for the cor- 
responding portion of equal mass flows, one obtains 


1 + (n/¢B)? 
For equal mass flows, Eq. 43 reduces to Eq. 14 when 6 = 1. When 
cos |{ 1 | 


Multilayer Rotating Wire Screen Heat 


4. Temperature Effectiveness of a 
I-xchanger 


Parallel Flows 


A procedure similar to that employed for equal mass flows yields 


=> Ye 


= m + 


= no + 


= 13 + 


For the special case of 9., = 162 = Me; = 1s, = Ns, Eq. 45 reduces to 


| 
| 
; 
4 
| 
1 f 
1 — 2 + 
1; 1 — 22 (5 +5 
(45) 
£13 
= + 1—2 = +> 
= Nsn Nn-1 9 9 
= 
Where 
W. 2-5 
- 


| 
= n2+7.|1 


= t+ 


= Ya-1 + 


b. Counter Flows 


A procedure similar to that employed for equal mass flows yields 


m1 = Ne 


+ 83 


= ). 


For the special case of 7., = 7.. = 


(75) 


= Ns 


= m1 


nenif 


1 nsnaf 


-2n1(5+4)] 


1—m 


1 Nef ) ( 


— %2 


; = n, Eq. 47 reduces to 


nensf 
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DESIGN PROCEDURES AND SAMPLE CALCULATIONS 


Equations developed in the Analysis Section may be readily plotted 
as design charts to establish design criteria. This was done (see Figs. 
8-11) for equal mass flows. For unequal mass flows, similar plots may 
be readily prepared, using the equations for unequal mass flows. These 
plots are not shown here, but are left to the individual designer to 
prepare to suit his own applications. For purpose of illustration, the 
general design procedure is given below: 


Gee Ta Pu Pa Wu Con Coe 
Required: To design a heat exchanger to attain an effectiveness 7. 
Procedures: 

1. Arbitrarily select the following quantities: m, a, n,, d, s, counter 
flow or parallel flow, and wire material which determine p,,C),,. 

2. Read n, in Fig. 8 or 9 from known values of n and ”,. 

3. Select ¢ and @ in Fig. 10 from known value of n,. (It is always 
desirable to have 6 as close to maximum as possible and ¢ as 
close to minimum as possible.) 

4. Calculate ( : s-) from Eq. 15 and read (“= ) in Fig. 11 
Pp 


ma 


from known value of (“- ). 


| 
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5. Calculate G from Eq. 13, N from Eq. 16. 

6. Calculate free flow area A; for either hot or cold gas from 
= —- = —~, and calculate the outside diameter of the wire 
Gi, 
screen rotor, using 


2 f 
0 


NOTE: (2) may be arbitrarily selected, but is usually from 0.1 
0 


to 0.3. 
Calculate the total wire screen core weight W’, using 


d \? 
W. 4 (D5) (>) 2s pwn. 


Calculate the core thickness of the wrire screen layers /, using 


| 


l= (n,-—1)p+d (in inches). 


9. Calculate the core volume using 


Core Volume (Po) (>: 


205 
+ 
| | | | 
oh 
5 - t—-+—ft +—— + 005 + 
| | | 
| | 
a | | 
| 
j 
| 
g.05 
| | 
02 + + + + + 4 + + + + 
| } | | 
01 
2000 3000 5000 19000 
R 
5 
LE 
: 
8. 


206 W. H. T. Lou (J. F. I. 


For a given rotating wire screen heat exchanger, the performance 
can also be calculated from the charts in a similar, but reverse, manner. 
In other words, knowing the values of V., N, etc., the value of 7 may 
readily be calculated for the given heat exchanger under the given 
operating conditions. 

SYMBOLS 


= free flow area in sq. ft. 
= specific heat of gas in BTU/Ib. °R. 
= specific heat of wire material in BTU/Ib. °R. 


= Constant = ( 


= inside diameter of rotating wire screen rotor in ft. 
= outside diameter of rotating wire screen rotor in ft. 
= wire diameter in ft. or inches 
We 
G ~~ = mass flow velocity per unit area in Ib./sec. ft.2 = pV 
Gmazr = Maximum mass flow velocity per unit area in Ibs./sec. ft.2 
g = acceleration of gravity in ft./sec.? 
h = heat transfer film coefficient between gas and wire in BTU/sec. ft.2 °R. 
k = thermal conductivity of the wire in BTU/sec. ft.? °R./ft. 
l = thickness of the rotating wire screen rotor in ft. or inches (thickness of multi-layer 


screen part only) 
M = Mach number = ( =) 
VygRt. 
m = total number of hot and cold sections in duct (see Fig. 12) 


a 


(mM = 4 IN THIS SKETCH) 


Fic. 12. 


revolutions of rotating wire screen rotor per second 
integral number = nt* harmonic in a Fourier Series 
number of screens in parallel in a rotating wire screen rotor 
total pressure in lb./in.? or Ib. /ft.2 
Prandtl number 

= gas constant in ft./°R. = 53.3 (for air) 


4 
: 
: 
: 
: 
x 
Ay 
Cp 
C 
Pw 
4h 
: 
; 
|_| 
N 
n 
if 
P 
P. 
r 
a 
By 
5 
t 


Sept., 1961.] 


General Subscripts 


Rotary Matrix WIRE SCREEN EXCHANGER 


Transverse pitch between wires in ft. or inches 


c indicates cold gas 

h indicates hot gas 

2 indicates initial 

n indicates n** harmonic in a Fourier Series 

o indicates steady state time wire element entering hot section 

s indicates single wire screen 

w indicates wire or wire screen 

x or bx indicates steady state time wire element leaving hot section 


T = total temperature in °R. 
T,;, = total temperature of hot gas leaving the wire in °R. 
T» = instantaneous wire element temperature in °R. 
AT. = wire temperature rise in °R. during passage of the wire through the hot section 
= Tor — 
t = time in seconds 
V, = axial velocity of gas in ft./sec. 
W = weight flow of gas in lb./sec. 
W. = total wire weight in lbs. 
a = angle in degrees or radians which prevents mixing of hot and cold gases during passage 
through rotating wire screen rotor 
_1/ eC, 1 Va 1 — ma/2r 
= 2( (aa) (na) 
Y = ratio of specific heats of gases 
h d/s 
n = over-all temperature effectiveness of rotating wire screen heat exchanger 
m. | = temperature effectiveness of an individual wire screen 
™1, 72, M3, *** Mm = temperature effectiveness of one, two, three, four . . . m wire screens in 
parallel 
Ney, Ney *** Neon = temperature effectiveness of first, second, third, fourth . . . . n*® indi- 
vidual wire screen 
0 = angle in degrees or radians in hot or cold section of the duct 
p = density of gas in lb./ft.3 
d d\? 
= screen solidity = 2(£) (£) 
w = equivalent angular frequency in radians per second = ~ ee 
1 — ma/2r 
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THE MINIMAL TIME REGULATOR PROBLEM FOR LINEAR 
SAMPLED-DATA SYSTEMS: GENERAL THEORY* 


BY 


C. A. DESOER' AND J. WING! 


ABSTRACT 


This paper considers a single-input, linear, time-invariant, sampled-data system 
described by its state transition equation x,,, = Ax, + f%,:a, where the control fx,1 
is subject to the admissibility condition | fx,;| < 1. Necessary and sufficient condi- 
tions for controllability with admissible controls are derived. For any system 
satisfying these conditions, the following problem is solved: determine a scalar valued 
function f(x) such that if, at each sampling instant, f(x) is used as a control when the 
system is at state x at that sampling instant, then, given any initial state Xo, the 
system will be brought to the origin in the minimum number of sampling periods. In 
other words, f(x) constitutes an optimal strategy for the minimal time regulator 
problem. The function f(x) is describable as follows: (1) in the state space, a hyper- 
surface @ is constructed; (2) the line through x, parallel to r, = — A™'a intersects @ 
at point c, thus x = c + Ar); (3) f(x) = sat A. 


I. INTRODUCTION 


In recent years considerable attention has been given to the problem 


of optimal time control for linear continuous systems (1-9).2. The most 
important aspect of the problem is that the control signals are restricted 
to belong to a closed and bounded set. In the case of a single control 
signal f, the restriction is that | f(t)| < y, where y is a prescribed num- 
ber. The fact that f is bounded plays a pre-eminent role in the design 
of single-input optimal systems. ‘The analysis of optimal control of 
linear discrete systems has been the object of much less attention 
(10, 11, 12). For discrete systems the additional restriction on control 
signals is that the values of the signals can change values only at fixed 
instants of time. The present paper solves completely the single-input 
minimal time regulator problem which was considered in special cases 
in two previous publications: in (11) a second-order case was considered 
and in (12), the m‘»-order case was restricted to transfer functions that 
have real, distinct, non-positive poles. It is complete in the sense that 
no restrictions are imposed on the location or multiplicity of the poles 
of the transfer function other than those necessary and sufficient condi- 
tions for the system to be controllable (13). 

* This research was jointly sponsored by the Electronics Research Directorate of the AF 
Cambridge Research Center, Air Research and Development Command, under Contract 
AF 19(604)-5460 and the Air Force Office of Scientific Research of the Air Research and 
Development Command, under Contract AF 18(1600)-1521. 

1 Department of Electrical Engineering, University of California, Berkeley, Calif. 

* The boldface numbers in parentheses refer to the references appended to this paper. 
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This paper considers a sampled-data control system described by its 
state transition equation X,,, = Ax, + f:,:a where A and a are con- 
stant and X1 matrices, respectively. Such an equation 
describes, for example, a system consisting of a sampler followed by a 
hold circuit whose output controls a black box with a transfer function 
G(s) which is a constant divided by a polynomial in s. A constant 
coefficient state transition equation such as the one above may describe 
much more general systems, for example, under certain conditions, 
delay lines or transportations lags are allowable. 

The minimal time regulator problem may be stated as follows: 
given any arbitrary initial state, it is required to bring the system to 
equilibrium in the minimum time with a control function f constrained 
by the condition |f| <1. It is furthermore required to find the 
optimal strategy that will generate such an optimal control for all 
possible initial states. 

The problem is precisely stated in Section II. In Section III, the 
general necessary and sufficient conditions for controllability with ad- 
missible controls are stated. Sections IV, V, and VI define and estab- 
lish properties of some sets Ry’ required to arrive at the proposed 
optimal strategy which is described in Section VII. Finally, in Section 
VIII, it is shown that, as T — 0, the proposed optimal strategy be- 
comes identical with the switching surface technique for generating the 
optimal control in the continuous case. In this paper we concentrate 
on establishing carefully the general theoretical framework. Since, in 
the general case, the optimal strategy will necessitate measurements 
analogous to those required in the special cases treated in (11 and 12), 
no description of possible implementation of the optimal strategy is 
included in the present paper because such descriptions are found in 
(11 and 12). 


Il. STATEMENT OF THE PROBLEM 


We consider a single-input, linear, time-invariant, sampled-data 
system. Let x, (k = 0, 1, 2,---) be the state of the system at time 
kT, where T is the sampling period. The input which will, henceforth, 
be referred to as the control is required to be constant during each 
sampling period. The control has the value f; for (k — 1) T < t < kT, 
k = 1,2,---. The system's state transition equation is 


k=0,1,2,-:- (1) 


where Ais an 7 X constant matrix, a is a constant column vector and 
x, is an n-rowed column vector. Since in many applications such a 
system results from the sampling of a continuous system, the matrix 
A = eB’ where B is the matrix appearing in the differential equation 
describing the continuous system. Since for any constant matrix B, 


ay 
| 
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det (eB?) ~ 0, we assume 
det A # 0. (2) 


In order to represent the common phenomenon of saturation, we as- 
sume throughout that 


lfe| <1 for all k. (3) 


All controls satisfying (3) are called admissible controls. The problem 
we propose to solve is referred to as the minimum time regulator problem 
and is the following : given an arbitrary initial state Xo, find an admissible 
control (specified by fi, fs, ---) which will bring the system to equi- 
librium (that is, the state x = 0) in the minimum number of sampling 
periods. Any such admissible control is called an optimal control. 

It is also required to find for every possible x» a method for generat- 
ing an optimal control. More precisely, it is required to find a scalar 
valued function of the state x, f(x), such that: (a) | f(x)| <1 for all 
states x; (b) if at each sampling instant, f(x) is used as a control when the 
system is at state x at that sampling instant, then given any initial 
state Xo, the system will be brought to the origin in the minimum number 
of sampling periods. The function f(x) is called an optimal strategy 
because it constitutes a rule for generating optimal controls irrespective 
of what the initial state x» might be. 


Ill. PRELIMINARY ANALYSIS AND CONTROLLABILITY 


Since no assumptions have been made concerning A and a (outside 
the fact that they are constant and (2)), we must first determine the 
conditions under which the problem has a solution. Suppose the con- 
trol fi, fo, ---, fx were applied, then 


x, = Ax, +/fia 
x, = A’x, + fiAa + fra 


Xx = Ax, + fiA*—a + + fa. 


Now let 
r, = — Aa (4) 


and we have thus established the fact that the system can be brought to 
the origin in N sampling periods (that is, xy = 0) if and only if 


(5) 


Let us first consider the problem when the control signal f; is not 
required to be admissible, that, is, it is not subject to inequality (3). 
We introduce now a concept due to Kalman (13). 
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Definition: The system (1) is controllable if an input can be found with 
the following property: given any initial state xo, the system can be 
brought to the origin 0 in a finite number of sampling periods. 


Theorem I (13). The system (1) is controllable if and only if a, Aa, 
A’a, ---, A”~'a are linearly independent. 


Proof: From (5), Xo can be brought to the origin in 7 sampling periods 
if and only if x» can be written as 


Xo =fiti t+ foe t+ 


Therefore, for any Xo, a control { f1, fo, ---, fn} exists with that property 
if and only if {r:, m2, ---, r,.} constitute a basis for the state space; in 
other words if and only if A~'a, A~*’a, ---, A~"a are linearly independent. 
Since det A # 0, this is equivalent to a, Aa, ---, A"~'a being linearly 
independent. 

Finally, if the system cannot be brought to the origin in » sampling 
periods, this cannot be done in any finite number of sampling periods. 
This is a consequence of the Cayley Hamilton theorem (14) which 
states that A* is a linear combination of I, A, ---, A"~'; therefore A’a 
is a linear combination of a, Aa, ---, A™—'a._ By iteration, it follows that 
A‘a (k > n) has the same property. Consequently if the system can- 
not be brought to the origin from some X» in ” sampling periods, it is 
impossible to do so in any finite number of sampling periods. This 
completes the proof of Theorem 1. 


In the above result, the controls were not required to be admissible. 
We consider next the effects of this additional restriction. 


Definition: A system is said to be controllable by admissible controls if for 
any initial state Xo, there is a finite length admissible control fi, fs, «++, fx 
which brings Xp to the origin 0. 

Observe that the duration of the admissible control, N, depends on 
X» and it is not bounded by a duration M independent of Xp as it is in 
the case where the controls are not required to be admissible. 


Theorem 2: The system (1) is controllable by admissible controls, if and 
only if 


(a) a, A'a, A’a, ---, Aa are linearly independent, 
(b) the eigenvalues of A are in the closed unit circle, (that is, 
< 1). 


Proof: Let us prove first, necessity. We assume that the system is 
controllable by admissible controls. Hence, it is controllable and (a) 


2 
: 
| 
3 
A, 
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holds by Theorem 1. Suppose now that (b) does not hold, that is, a 
particular eigenvalue is outside the circle, say |\,| > 1. By assump- 


tion, with N sufficiently large and the |f;| <1, the sum — > fA-‘a 

can represent any point in the state space. Consider the fundamental 

formula® for A~‘ and the resulting coefficient of Z,oa of the above sum; 

it turns out to be 


N 
firs 
i=l 


Now, however large N is, 


Therefore, no admissible control exists which will make the coefficient 
of Z,oa larger than 1//A,| — 1. Hence, for example, if x» = &Z,oa with 
€>1/|dAx| — 1, Xo is not controllable. This contradicts the assump- 
tion, therefore, no eigenvalue may be outside the unit circle and (0) is 
established. 

Let us next prove sufficiency : we assume that (a) and (b) hold and 
establish controllability by admissible controls. Since (a) holds and 
since det A 0, all the sets of vectors (fpn41, *)Tpnin), for 
p = 0, 1, 2, ---, form a basis for the state space. Given any X» we 
define x”, for p = 0, 1, ---, as follows: expand X» with respect to the 
p* basis, thus 


n 


i=] 


For each p, divide each n;‘”) by max |n;‘”|; let the result be &;; by 
1<i<n 
construction |&;(")| < 1. The vector xo‘” is defined by 


(6) 


The vector Xo‘”? lies in the direction of X» and has positive (>0) length. 


Clearly Xo is controllable by admissible controls in at least Nn sampling 
N-1 


periods if has a smaller length than x» and > has a length 
p=1 

greater or equal to X»._ It remains to show that given any Xp» such a 

finite N exists. The only way in which no such finite N would exist 

is for as p— This can only happen if © as 

pb — « forsome 1 <i <n. Considering the expression for 7; given 


See Appendix I. 
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by Cramér’s rule in solving (6), this can only happen if 

det (fpn41) Tpnt2, °° 7 O as p- 
However, since = 


det (fpn41, °° = (det det (11, Fn). 


Now det A = JJ \,” where m, is the multiplicity of the eigenvalue ),. 
k 


=1 
From (b), |det A] <1 hence xo‘”) cannot go to zero as p — © and 
the existence of N is established. Q.E.D. 


Since the problem stated in Section II is meaningful only if the 
system is controllable by admissible controls, henceforth we make the 
following 


Assumption: The matrix A and the vector a are such that 
(a) a, Aa, A’a, ---, A”~'a are linearly independent (7) 
(b) the eigenvalues of A are in the closed unit circle (that is, 
< 1). 


A direct consequence of this assumption and of Eq. 4 is that: for any 
integer k, the m vectors fy, Te41, Teyn—1 are linearly independent. 


IV. THE SETS Rw’ AND Rw 
It is obvious that since we are interested in generating for each 


initial state an optimal control the sets Ry’, to be defined presently, 
will play a leading role. By definition, for any integer JN, 


Ry’ is the set of all initial states that can be brought to the origin (8) 
by an admissible control in N sampling periods or less. : 


It follows immediately that Ry’ defined by (8) is characterized by 


Property 1: (1) Ry’ is also the set of all states that can be written 
as a sum 
m0, + where |n;| < 1 for? N; symbolically, 


N 
Ry’ = {x|x = |ni| < 1 fori = 


=] 
(2) Ry’ is convex, that is, if x,;€ Ry’ and x. € Ry’, then x = Ax, 
+ (1 — A)x. € Ry’ for all AE [0, 1]. 


Suppose, for example, that m = 3 and that rm, fo, r;, and ry are as 
shown in Fig. 1, then R,’ is the segment A,A, and R,’ is the parallelo- 
gram B,B,B;B, shown in Fig. 2; R;’ is the parallelepiped C,C2, ---, Cs 


| 
| 
| 
| 
| 
| 
| 
| 
| 
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shown in Fig. 3; Ry’ is the prism D,Dz, ---, Die shown in Fig. 4. Con- 
sider another example: n = 2 and the eigenvalues of A are exp(+ ja) 
where a = 53°; the vectors fr, f2, «++, fs are shown in Fig. 5; the sets 
R,’, R:', ---, Rs’ are shown in Fig. 6. 


A2 
3 B, 
hy 
B, 
| A, 
Fic. 1. The vectors r, fe, f3, and ry, of Fic. 2. The parallelogram Ry’. 
the example illustrated in Figs. 2-4. 
D, 
| 
| Cs 
D, | 
C, | 
| Cc 
a] 
| 
4 | \ 
Jj. D, 8 
Cy — — Ce 
4 
/ “Cs 
Cy Cs De 
Fic, 3. The parallelepiped R;’. Fic. 4. The polyhedron R,’. 


For future use we have to introduce the idea of a convex hull (18). 
By definition, the convex hull of a finite set of points S = {y:, ye, ° *, Y>} 


is the smallest convex set containing S; or equivalently, y = ¥ \Xy; 


is a point of the convex hull of S if and only if (a) A; > 0 and 


3 | 
‘ 
hy 
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The vectors fr, f2, fs for the Fic. 6. The regions Rs’. 
example illustrated in Fig. 6. 


; = 1. Physically the convex hull of S, denoted by [S], is 


the set of all points y which are centers of gravity of the points y; of S 


when y; is assigned a mass )j. 
N 


Let, now, Vy be the set of 2% points given by >} et; where the e; 


can take either the value —1 or 1. We then have 


Property 2:4 Ry’ is the convex hull of Vy; symbolically, Ry’ = [Vw]. 
This property is interesting because it specifies that Ry’ is a polyhedron 
and that its vertices are points of the set Vy, although it does not 
specify which of these are vertices and which are interior points of Ry’. 
Now for control purposes, in addition to Ry’ we are interested in Ry, 
defined as follows: 

Ry its the set of initial states that can be brought to the origin by an 
admissible control in N sampling periods and no less. Ro is the origin 
{0} itself. 


N 
Clearly from the definitions, Ry’ = U R:, that is, the union of 
=0 


Ro, Ri, Ro, «++, Rv. Also Ry = Ry’ — R’yw-i, that is, Ry is obtained 
from Ry’ by deleting from it all points of R’y_;. A useful relation 
between Ry and R’y-_, is given by 


Property 3:4 lf x» € Ry, then Xo = p + nyFw where p is a boundary 


point of R’y_; and |ny| <1. 
This property calls for two remarks: 


4 The proof of this property is in Appendix II. 
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(1) The converse to property 3 is not necessarily true, that is, 
if’ and |ny| <1, x = p+ nary may not bein Ry. This 
may be the case either because the vector nyfty points (from p) into 
R’ y_, or because ry lies in the face of R’y_, that contains p and ny is 
small enough so that p + nary is also in that face. 

(2) There are some cases for which the representation x») = p 
+ nyTy is not unique. 


V. CONSTRUCTION OF R’y,:1 FROM Ry’ 


From properties 1 and 2 it follows that, for any integer NV, Ry’ 
convex polyhedron. In particular: 


R,’ is the line segment, centered on the origin, joining r, to —r,, that is, 
= {x|x = fin, | fil 

R,’ is the parallelogram, centered on the origin, whose edges are parallel 
to rm, and fp; vertices are —f, — + fs, 1 + £2, — 
Symbolically, = {x|x = fir, + forts, | fi] <1, | feo] < 1}. 


R,’ is the hyperparallelepiped, centered on the origin, whose edges are 


n 


parallel to r,, ---,r, and whose 2” vertices are of the form © er, 
1 


where the sequences (e;) constitute the 2" possible sequences of ones 
and minus ones. 


For N > n, Ry’ is obtained inductively. Suppose Ry’ is known, let us 
construct R’y,;. Ry’, being a convex polyhedron, is completely de- 
fined by its vertices (see property 2). Let us classify the vertices of 
Ry’ in two classes: in the following 6 is an arbitrarily small positive 
number. The class Vy* contains all the vertices p; of Ry’, such that 
Pi + 6ry,,€ Ry’. The class Vy~ contains all the vertices q; of Ry’, 
such that q; — 6ry,,€ Ry’. These two classes usually have common 
vertices, see Fig. 4, the vertices Cs, Cy, Cs, Cs. It follows from the 
fact that, for all N, Ry’ is a convex polyhedron and from properties 2 
and 3 that we have 


The vertices of R’y., are obtained from those of Ry’ 


Property 4 


as follows: 


(a) to each vertex p; of corresponds the vertex of R’y 41: pi + 
(b) to each vertex q; of V y~ corresponds the vertex of R’y41:qi —Tw4t 


This property follows directly from properties 1, 2, and 3. Let us, 
now, comment on the geometrical aspects of the re ation between 


5 We use the eaten’ notation of dRy’ to indicate the naiany of Ry’, that is, the set 
of all points x, every neighborhood of which contains at least one point of Ry’ and one point 
of the complement of Ry’. 
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R’y,, and Ry’. With respect to the direction ry,, (@) Rwy’ has a 
“top,’’ namely, those faces that have all their vertices in Vy*t, (6) Ry’ 
has a ‘‘bottom,”’ namely, those faces that have all their vertices in V y-. 
OR’ the boundary of R’y 41, is obtained by raising the ‘‘top”’ of Ry’ 
by and lowering the “bottom” of Ry’ by —ry4:. 

As an illustration consider the case shown in Figs. 1-4. By con- 
sidering Figs. 1, 3, and 4, it follows that V;+ consists of Ci, Ce, Ce, Cs, Cs, 
and C4; V3~ consists of Cs, Cs, Cz, Cs, Ca, C3. Note that the vertices 
Co, Cs, Ca, Cs of Rs’ belong to V;+ and V;~. This example is a particular 
case since ry is parallel to the faces Ci, Cs, Cs, Ca and Cs, Ce, Cz, Cs of 
R,'. Asaresult the face D,, Ds, D3, Ds, Ds, Ds of Ry’ contains the face 
Ce, Gs, of R,’. 

It is intuitively obvious that 


Property 5:6 For each N, Ry’ has edges parallel to 1), fe, «++, Py. 


VI. DEFINITION AND CONSTRUCTION OF THE CRITICAL SURFACE 


The critical surface @ is an (x — 1)-dimensional hypersurface which 
will play a crucial role in obtaining the optimal control. We shall 
construct it step by step in a process to be described presently. 

Let us consider a sequence of sets of points: p; contains only the 
origin. ps is the segment joining to re, that is, ps = {x|xX = 
<1} +++ px is an — 1)-dimensional parallelepiped defined by 


n 


Let us note that since the vectors are linearly inde- 
pendent, the sets p:, ps, +++, pn lie in the (7 — 1)-dimensional subspace 
spanned by fe, ++, 

For any N > n, we define py,; in terms of py and Ry’ inductively. 
Consider the intersection of py with the boundary of Ry’, that is, the 
set py (1) ORy’. Since Ry’ is a convex polyhedron this intersection 
consists of sets of points that lie on faces of Ry’. Let these faces be 
denoted by ®y.1, ®y,2, To construct py,; we consider each face 
@y,; and (1) if ry4; is not parallel to any vector lying in the face ®y |, 
we construct the set Vy,;,; defined by 


= = y + Ary,, where y € Py,; pw and X belong 
either to (0, 1 ]or [—1, 0) the choice being made sothatx€ Rwy’} (9) 


(2) if ry. is parallel to some vector of the face ®y,;, then Wy41,; 
= py) ®y,;. Note that in case 2, Vyii1,; C py. By definition, 


= U V (10) 


6 This property is proved in Appendix II. 
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This construction is illustrated by considering again the example of 
Figs. 1-4. Figure 7 shows the parallelepiped R;’ and ps which is the 
parallelogram Fi, F2, Fs, Fs. The faces P32, P33, and are indi- 
cated. Consider now Fig. 8. Since ps; is the line segment Fi, Fs 


w 


Co Ge Ce 


c 
8 


Fic. 7. The parallelepiped C;, ---, Csis R;’ Fic. 8. Illustration of the construction of p, 
of the example of Figs. 1-4. from p3. 


and since r, is not parallel to any vector of the face ®;,, the parallelo- 
gram Wy, is shown as F,, Fy, G:, Gs. Similarly for the face 3; the 
parallelogram V4; is shown as F2, F;, G2, G3. Now, since r, is parallel 
to a vector of the top and bottom faces of R;’, Was is the segment 
F,, F, and V4, is the segment F;, Fy. p3 consists of two parallelograms 
Fs, Gi, Ga; Fe, Fs, Go, G3 and two line segments F;, and F;, 
Reference to Fig. 6 of (12) shows that p, consists of the parallelograms 
pest, Psat, Pos, ps Of the parallelograms post, Psst, Past, 
and 

Before defining the critical surface, let us add a few remarks 
on the sets py. (a) Since py is an (n — 1)-dimensional set, the points 
y © py 1) ORy’ form an (n — 2)-dimensional set. Therefore, in case 1, 
the set Wy41,;, which is obtained by sweeping all points y € py () ®n,; 
in an outward direction along the vector ry4,, isan (m — 1)-dimensional 
polygon. In case 2, Wy,,,; is an (m — 2)-dimensional polygon. In 
either case, the polygons Vy,,,; are connected to py. In either case, 
also, the Vy,,,; start on a face of Ry’ and terminate on a face of R’ w41. 
pai, the union of all the Vy,,,;, is a set of polygons that join faces of 
Ry’ to faces of R’y4:. (b) Observe that py C Ry. 


Definition of the Critical Surface C: C = U p, 


‘ 
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A useful property for the derivation of the optimal control is 


Property 6:7 For each integer N, the critical surface @ intersects the 
boundary of Ry’ only on faces containing edges parallel to r;. 
This allows us to say that py4, is the union of polygons starting on 
faces of Ry’ containing edges parallel to r,; and terminating on faces of 
R’ v4; containing edges parallel to r;. See also Fig. 8. 


As a final step we state 


Property 7: Any straight line parallel to r, intersects the critical 
surface at one and only one point. 


Let us stress that property 7 states that for any such straight line 
an intersection point exists and that it is unique. 

Referring again to the example of Figs. 1-4, we note that the portion 
of the critical surface contained in R,’ is the two-dimensional surface 
Go, Fs, Fi, Gi, Gs, Fs, F3, Gs shown in Fig. 8. In the example of Figs. 
5 and 6, the critical surface is one-dimensional, hence is a polygonal 
line. The portion of it that lies in Rs’ is shown in Fig. 9. It is easy to 
observe that in both cases properties 6 and 7 are satisfied. 


Fic. 9. This illustration refers to the example already considered in Figs. 5 and 6: the 
heavily drawn polygonal line is the portion of the critical surface @ (which is one-dimensional 
in this case) which lies in Rg’. 


7 The proof is in Appendix II. 
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VU. THE PROPOSED OPTIMAL STRATEGY 


We can now reformulate the definition of optimal control (given in 
Section II) in terms of the concepts Ry’ and Ry. Suppose the initial 
state Xo © Ry, then the admissible control f; is optimal if and only if the 
state x, R’y_1, where x, = Ax, + fia. Observe that x» Ry and 
x, € R’y_, imply that x; Ry_;. An optimal strategy is a rule which for 
any state x gives an optimal control for that x; it is a function whose 
domain is the state space and whose range is the set of admissible con- 
trols. Any computer that generates at each sampling instant an 
optimal control is simply a realization of an optimal strategy. Con- 
versely the first step in the design of such a computer is to find an 
optimal strategy. 

First let us introduce a 


Definition of the Critical Slabs § 


The set ¢ consists of all points x such that x = c + Ar, where Cc isa 
point of @ and |A| <1. Symbolically, 


[al <1}. (11) 


Note that for a given x, ¢ and \ are unique, by property 7. 
We shall derive the optimal strategy in successive steps: 


(1) If the initial state x» is on the critical surface @, then f; = 0 
is an optimal control for the first sampling period. This is immediately 
established : since X» € @ and there is an N such that xX) € Ry, we have 

N 
Ry @, hence xX» has a representation of the form = na; 

where |;| <1 for 7 = 2, 3, ---, N. Suppose we apply the control 
f, = 0 during the first sampling period. At the end of it, the system is 
in state x, and by (1) and (4), 


V-—1 
x, = Ax, = 
i=1 


Hence, by property 1, x, © R’y_; and f; = 0 is optimal. 
(2) Suppose the initial state xX» is in the critical slab ¢, then f; = A 
is an optimal control for the first sampling period, where \ is defined by 


X» =c+OAr, with cE e@. (12) 
Let N be such that x»€ Ry. Since py, hence ¢, intersects 0Ry’ on 
faces containing edges parallel to r; and since Ry’ is convex, the point ¢c 
defined by (12) is also in Ry’, hence 


N 
Xo Ar, + "iT; 


9 


i= 
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where |n;| <1 = 2, 3,---,N. Sincexo€ ¢, |A| < landf; =A 
is an admissible control. Again using (1) and (4) we have 


N-1 


x, = Ax» + Ab = 


where < 1 ford = 1,2,---,N—1. Hencex,€ R’y_,andf; => 
is optimal. 

(3) Suppose X» does not belong to the critical slab ¢, then f; = sgn 
is an optimal control for the first sampling period, where \ is defined 
by (12). 

Let N be such that x»€ Ry. Since xo € ¢ and x, € Ry’ thec and \ 
defined by (12) are such that c€ Ry’ and |\| > 1. Therefore the 
lemma of Appendix III applies and 


N 
Xo = (sgnaA)ri + 


i=2 


with |y;| < 1,7 = 2,---, N. Therefore, as above, f; = sgn \ is both 
admissible and optimal. 
These three results allow us to formulate the 


Proposed Optimal Stragegy 


Suppose that at the present sampling instant the system is in state 
X». Compute \ such that 


Xo 


wherec € @. Then/f,; = sat \isan optimal control for the next sampling 
period.’ The function sat ) is the function f(x) mentioned in Section IT: 
it is the optimal strategy. 

Let us note that the above stated optimal strategy is not the only 
optimal strategy. For an illustration of the nonuniqueness of the 
optimal strategy see (11), Fig. 10 in particular.2. The proposed strategy 
seems to be the simplest and also the one that will give naturally the 
optimal strategy for the continuous case as T — 0. 

Let us note also that the proposed optimal strategy requires only 
the simulation of the critical surface @ and the determination of \, 
the distance along the direction r,, of the present state x to the critical 
surface @. For special cases, examples are to be found in (11 and 12). 


8sat x = x for |x| < 1 and sat x = sgn x elsewhere. 
® Another optimal strategy would then be as follows: at xveRy, take f; = \’, where X’ is 
the A smallest in absolute value such that Axo + fiaeR’y_1. 
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VII. RELATION BETWEEN THE CONTINUOUS CASE AND THE DISCRETE CASE 


In this section we propose to discuss heuristically the relation be- 
tween the optimal control of a continuous system and that of the corre- 
sponding sampled system. More specifically, we consider a fixed con- 
tinuous system 


x= Bx + f(i)b (13) 


where B and b are, respectively, m X » and n X 1 constant matrices 
and |f| <1. Associated to the continuous system (13), we consider 
the discrete system obtained by sampling the control f(t) and keeping 
it constant during the sampling period: thus, for the sampled system 
the control is f, = f[(k — 1)7] for the interval (k — 1)T <t < kT. 


The resulting discrete system is characterized by 


Xigi = A(7)x, + fixsa(T) (14) 


where the dependency of A and a on 7 is explicitly indicated. The 
following relations will be useful. 


(15) 
r T? T3 
a= = Tb — 7 Bb + Bb — (16) 


It is known (7) that, for the continuous system, if b, Bb, ---, Bb are 
linearly independent and if the eigenvalues of B have negative real 
parts, then, for each initial state Xo, (1) there exists an admissible 
minimal time control that brings the system to 0, (2) the minimal time 
control is piecewise constant and (3) the set of all admissible controls 
that simultaneously satisfy the maximum principle and bring the 
system to 0 consists of only one element. 

In order to show how the control and the optimal strategy of the 
discrete system merge into that of the continuous system as 7 — 0, we 
consider a sequence of sampling periods: 7; = T/2!,/ = 1, 2, --- and 
T) is a fixed number. For each /, the proposed optimal control will 
bring the system to the origin starting from X»in a time 7(Xo, /). Clearly 
7(Xo,/ + 1) < r(Xo,/) because the optimal control for 7; is still an 
admissible control for = 7T,/2. Therefore the sequence of 7 (Xo, /) 
forms a monotonically decreasing sequence of positive numbers, hence 
as ~, 7(Xo,1l) 7(xo). The limit 7(xo) is equal to 17.(Xo), the 
minimal time required by the continuous system (13) to bring the 
system to 0 starting from x». This follows because, for sufficiently 
large /, a piecewise constant (over intervals of length 7;) approximation 
of the optimal control of (13) will be an admissible control of (14) and 
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the corresponding solution of (14) is an ¢e-approximation (16) to the 
solution of (13). Therefore, as ] — , this procedure will give a solu- 
tion of (14) which will require a time arbitrarily close to 7.(xo). There- 
fore 7(Xo) cannot be larger than 7,(Xo). Conversely, it is impossible to 
have 1(Xo) < 1.(Xo) because the optimal control of (14), for any /, is an 
admissible control of (13). Hence 7(xo) = 

Consider now what happens to the critical surface (which we shall 
denote by e(7;) to emphasize its dependence on 7;) as]— «. As 
1 — ow, the critical surface e(7)) becomes a smooth surface which we 
shall denote by e@(0). We justify heuristically this assertion by the 
following observations: (1) the length of the edges of the polygons that 
make up @(7,;) go to zero as! — », because a — 0 as seen from (16); 
(2) the angle between successive r,’s ~0 as ]— ©. This can be 
checked by showing that the ratio of — 1%)/||r.|| ~O asl— 
even if k is increased so that kT; remains constant. The only case 
where the normal to the tangent plane of @(0) will not be a continuous 
function of position is at points where, in order to construct py41, the 
sign of ry,; had to be reversed with respect to the sign of ry used to 
construct py at that point. 

Since r; ~ 0 as/] — ~, the limiting behavior of the optimal strategy 
is clear: the sign of the forcing function is reversed as soon as the state 
space point reaches @(0). In other words @(0) becomes a switching 
surface. Thus we see that, as the sampling period approaches 0, the 
proposed optimal strategy merges into the well-known switching sur- 
face technique which is used to generate the optimal control for con- 
tinuous systems. 
IX. CONCLUSION 


This paper presented a complete derivation of an optimal strategy 
for the minimal time regulator problem for discrete systems described 
by equations of the form (1) and for which assumptions 2 and 7 hold. 
It has been shown that unless these assumptions are satisfied, the prob- 
lem does not make sense. The derivation is general in the sense that 
only the assumptions 2 and 7 are imposed on the system (1). These 
assumptions restrict the eigenvalues of A to have their magnitude 
positive and not greater than unity. The procedures described in the 
paper are applicable to all such cases irrespective of whether the eigen- 
values are distinct or multiple, real or occur in complex conjugate pairs. 
The proposed optimal strategy defines for every point of the state space 
an optimal control, in other words, it defines a function f(x) for all 
states x such that f(x) is an optimal control. From this point of view 
it is closely related to the concept of switching surfaces used in the 
continuous case (see discussion in Section VIII) and does not require 
consideration of the usually exponentially increasing solutions of the 
adjoint system introduced in (6 and 7). 
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APPENDIX I: THE CANONICAL EXPRESSION OF A‘ 


The purpose of this appendix is to state briefly the fundamental formula (17). 
Let A be a matrix with ¥(A) as the minimal polynomial and \j, A», «++, A, as eigenvalues. 


Expand TEN) in partial fraction and collect terms 
Ckl 


s 
1 as 


where m;, is the multiplicity of A, as a root of ¥(A). Define the matrices 


= (A — (A) 


where 


+++, A, then 


Now it can be shown that if /(A) is analytic in an open set containing Aj, Ao, 


In particular, for any integer 7, 


s 
g 


Ata 
k=1 1-0 


(A*) 
| 


This is the fundamental formula for A~ used in the proof of Theorem 2. 


APPENDIX II: PROOF OF THE PROPERTIES 1-7 


The purpose of this appendix is to give proofs for the properties given in the body of 


the paper. 


Property 1 follows immediately from the definitions. 
Property 2: Ry’ is the convex hull of Vy; symbolically Ry’ = [Vw]. 

Proof: In order to show that Ry’ = [Vy] we first show that (1) [Vw] C Ry’ and (2) 
Ry’ C [Vw]. The desired equality follows then from the definition of the equality of two sets. 
(a) It is immediate that R,’ = [V,] and, that, for any N, [Vw] C Ry’. The latter 


follows from that any point of Vy, that is, 2 ¢; rj, isin Ry’ by property 1. 


(b) The property will be established if the following induction is proved: [Vw-1] D R’w-1 
N 


implies [Vy] D Ry’. Let x be any point of Ry’ thereforex = 2 nx; with |ni| <1. Now 


| 


> nit; € R’y-1, hence by the induction assumption 


mri = Aspe 
i=l k~1 
where pi, are the points of Vy_1, Ax > 0, Ax = 1. Observe that ry is 
k=1 


a point of Vy, hence 


N-1 oN 
. 

X= 2 mii + = & Ape + 
1 k=1 


i- 


oN-1 
(5+ yz Ax( Pi +r) +(5- — Ty). 


ll 


| | 
ite 
aks 
: 
i 


Sept., 1961.] MINIMAL TIME REGULATOR PROBLEM 225 


This last equality shows that x is a point of [Vy]. Hence, the induction is established, hence 
[Vw]D Ry’ for all N. This together with the result of (a) proves the theorem. 


Property 3: lf Xo Rwy, then = + p where pis a boundary point of R’y_, and |ny| < 1. 
Proof: Since Xo € Ry, Xo has a representation of the form 


N 
x= with <1. 


Since such a representation is not necessarily unique, consider all possible such representations 
and, in particular, the one that has the smallest 7y in absolute value. Let that representa- 


tion be 


N 
Xo = ni'Ti. 
i=l 


Clearly we have |n;’| < 1,7 = 1, 2, ---, N. Obviously, 7x’ # 0 for otherwise x» would be in 
R’y-1 which is impossible since it isin Ry. Also 


N-1 
1 


N-1 
The point 2 nit; is a boundary point of R’y_,: indeed, let 7y’’ be obtained by reducing 7y’ 


tl 


by ein absolute value. 


The point Xo — 2y’’ry is not in R’y_, by definition of 7y’ but is arbi- 
N-1 


trarily close to the point 2 it; of R’w-1. 
tl 


Therefore, the latter is a boundary point (15) of 


R'w-1. 
Property 4 is established in two steps: (a) R’n41 D the convex hull of the vertices p; + ry41 
and q; — rwy: and (0) this same convex hull includes R’y,;._ The proof is immediate if one 


uses properties 1, 2, and 3. 


Property 5: For every integer N, Ry’ has edges parallel to rm, fo, «++, ry. 


Proof: The property is true for N = 1. Let us establish the induction: suppose it is true 
for N, prove that it is true for N + 1. By property 4, each edge of Ry’ parallel to r; is trans- 
lated either by ry,; or —ry4; to become an edge of R’y41, therefore R’y,; has at least one edge 
parallel to r; for 7 = 1, 2, ---, N. If rv41 is parallel to one of these r,’s the property is estab- 
lished. But if ry4, is not parallel to any of the r;’s (¢ = 1, 2, ---, N) then Vy* and Vy~ have 
at least one point, say q, in common. Let by definition q€ Vy*( Vy-, then any 6 > 0 
however small, q — 6rv4i1 & Ry’ andq + éry4: & Ry’. Therefore to this q will correspond an 
edge of R’ v4; parallel to ry41, namely the edge that joins q + rvy1 toq —rvy:. This estab- 
lishes the induction and property 5 is established. 


Property 6: For any integer N, the critical surface @ intersects the boundary of Ry’ only on 
faces containing edges parallel to r,. Observe that the intersection of @ with dRy’ is equal to 
that of py with dRy’. For 1 < N <n, the property is easily verified since from the above 
description of the Py and property 1, we have 


Ry’ = {xlx=y+ n, [Al <1, for 1<5N<n. 


For N > n we use induction: we assume that for some N > n, Py intersects Rw’ on faces 
containing edges parallel to r, and we have to prove that this is still true for Py41. 

Let @y,; be a face of Rw’ containing points of py. By the induction assumption ®y,; 
has an edge parallel to r;. Now case 1: ry4: is not parallel to any vector of y,; (that is, 


ry,, cannot be expressed as a linear combination of the edges of the face @y,;). Then all points 
This new face 


of @y,; undergo a translation of either or to become a face of R’y 41. 


: 
: 
: 
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will therefore also have some edges parallel;to r;. Now, the points y € #y,; () Pw will undergo 
the same translation to form the polygon Wy4:,; and the intersection of Py4, with OR’y,, will 
include points on the new face obtained by translation from @y,;. Therefore py,; intersects 
OR’ yy: on faces that contain edges parallel to r; and the induction is established for case 1. 
In case 2, fy41 is parallel to some vector of the face @y,;. After translation by ry, in an out- 
ward direction the new face of R’y,4; includes @y,; and Py; intersects this new face since 
Wy41,i has only points on this face. The induction is thus also established in case 2. 


Property 7: Any straight line parallel to r; intersects the critical surface @ at one and only 
one point. 
Proof: Let 1 be the arbitrary straight line, more precisely, let 
1 = w€(—%, ©), xo arbitrary but fixed}. 
(a) Let us first show that / cannot intersect Py at two distinct points a and b. 
For 1 < N <2, this is impossible because py is a set in a subspace spanned by 


(fo, and, by assumption, r, is linearly independent of the set of vectors r3, In). 
Therefore the property follows for N < n. 

Let us now use induction to prove it for N > m. More specifically we assume that the 
property holds for some N > n and wish to prove that it holds for N +1. We establish the 
induction by contradiction: suppose / intersects Pv, at two distinct points a and b. Clearly 
if rv4: is parallel to r;, reference to the construction of Py4; from Py shows that Pv4; C Py in 
this case; hence we have a direct contradiction with the induction assumption. It remains, 
therefore, to establish the induction for the case where ry; is not parallel to r;. With an 
obvious notation, 

a — b = (ua — 
Since a and b are in py41, by construction of py41, there are two points y, and yy such that 


a= + Ya E PN a) ORy’, |ra| < 1 
b = yo + ys € Pn) ORy’, <1 
Let us project these two relations along a direction parallel to r; onto an (m — 1)-dimen- 
sional subspace. Let ‘‘~’’ indicate the projected quantities. We have 
a=b 
Ja =4— 
=a— 
and fya1 #0. Since ya and ys C pw) ORy’, they lie on faces of Ry’ which contain edges 


parallel to r:, therefore f4 and §, are on the boundary of the projection Ry’. Equations A-1 
mean that §. and f are on a line originating from 4 and whose direction is that of fy4;. Since 


(A-1) 


Ry’ is convex, this line intersects dRy’ at two points p and q, where is the intersection point 
closest to 4. Since to obtain a and b from y, and ys one had to move these points along ry41 
in a direction outward with respect to Ry’, it follows that f. and §, must coincide with f. 
Consequently ye — Ys = vf: for some ». Note that » ¥ 0, for otherwise since ry,4; is not 
parallel to m1, we would also have a = b. Thus 


Ya — Yo = 77), vy #0 


which implies that there is a straight line parallel tor; (namely, x = ys + Ari, AG (— ~, @)) 
which intersects Py at two distinct points y, and ys. This contradicts the induction assump- 
tion. Therefore the induction is proved and the statement (a) proven for all N. 

(b) Let us next show that / cannot intersect € at a point a of py and at a point b of py, 


where b & py and N’ # N. 
Let us project orthogonally all the Rx’’s, px’s and / on a two-dimensional plane II perpen- 


: 
: 
: 
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dicular to r;. Indicate the projected sets by “~.”” Now isa single point Xo, Ry’ is a convex 
polygon of the plane I. For any N, p; C Ry’ for i < N implies that p; C R,’ for i < N. 
Observe that for any N, Pw is the union of polygons Wy,; which start on a face of R’y_, and end 
on a face of Ry’; each of these faces has some edges parallel to r;, hence project on II as a 
point or a line segment ; each of these points or line segments is on the boundary of Ry’. There- 


fore Py is outside the interior of R’y_;. Since By C Ry’, then By C Ry’ — (R’y-1).% Fur- 
thermore, reference to the method of construction of Ry’ from R’y_, and that of Vy,; shows 


that we actually have By = Ry’ — (R’y-1), for all N. This fact together with the general 
inclusion C R’y_; C Rw’ C C shows that if N # N’, by and by, cannot have 
interior points in common (although if |N — N’| = 1, they may have boundary points in 
common). The statement (6) follows from this last fact. 
(c) Any line / intersects @: This follows from the above reasoning because the assumed 


controllability by admissible controls of the system implies that U_ R,’ fills the whole state 
k=1 


space. Therefore U R,’ fills the whole plane Il. Therefore there isan N such that ¥ € Ry’ 


N 
and since Ry’ = U By, there isa k < N such that Xo € Ax. 
k=1 


APPENDIX III: LEMMA 


Lemma: Suppose fr, f2, «++, tw are linearly dependent. Suppose y has the following repre- 


sentations 
N N 
= + 2 at; = Bit, + 2 (A-2) 
where 
a;’>1 and lai] <1 fori = 2,3,---,N (A-3) 
and 
<1 for = 1, 2, ---, N. (A-4) 


Then y has a representation of the form 


N 
yernt+ 2 vr with < 1, 4= 2,3, N. (A-5) 


Proof: lf 8; = 1 the theorem is proved. Let us then consider the case where #; < 1. 
From (A-2) it follows that 


N 
(Bi — + = (8; — = 0. (A-6) 
Pick the number & as 
a’ —1 
— Br 


Observe that 0 < k < 1 because of (A-3) and (A-4). Now multiply (A-6) by & and add it 


to (A-2), 


N 
y = + — ay’) + 2 + k(Bi — ax) 


Hence 


N 
y=r+ iti 


10 The superscript ° over the symbol of a set indicates the interior (15) of that set. 


: 

; 

; 
: 
: 
: 
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: 
: 
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' 


228 C. A. Desorr AND J. WING 


where 
yi =ai + R(Bi — ai) = (1 — + 
Since, for i = 2, 3, ---, N, both a; and 8; belong to [—1, 1] and 0 <k <1 then 
lvi| <1 
and representation (A-5) is established. 


Corollary: Under the same assumptions and notations but for the first inequality (A-3) re- 
placed by a,’ < — 1, the conclusion becomes 


N 
with lys| <1 for i = 2, 3, N. 


REFERENCES 


D. C. McDonatp, “Nonlinear Techniques for Improving Servo Performance,” Proc. 
Natl. Electronics Conf., Vol. 6, pp. 400-421 (1950). 
A. M. Hopktrn, ‘‘A Phase Plane Approach to the Design of Saturating Servomechanism,” 
Trans. AIEE, Vol. 70, Part I, pp. 631-639 (1951). 
D. BusHaw, ‘Optimal Discontinuous Forcing Terms,”’ in Contribution to the Theory of 
Nonlinear Oscillations, Vol. 1V, edited by S. Lefshetz, Princeton, Princeton University 
Press, 1958, pp. 29-58. 
I. BOGNER AND L. F. Kazpa, “An Investigation of the Optimum Switching Criteria for 
Higher Order Servomechanisms,’’ Trans. AEE, Vol. 73, Part II, pp. 118-127 (1954). 
R. BELLMAN, I. GLICKSBERG AND O. Gross, “On the Bang Bang Control Problem,” 
Q. Appl. Math., Vol. 14, pp. 11-18 (1956). 
C. A. Desorr, ‘“‘The Bang Bang Control Problem Treated by Variational Techniques,” 
Information and Control, Vol. 2, pp. 333-348 (1959). 
G. BoLttyansk1, R. V. GAMKRELIDZE AND L. S. PoNTRYAGIN, Theory of Optimal 
Processes I,’’ Izvestia Akad. Nauk SSSR Seriya Matem., Vol. 24, pp. 3-42 (1960). 
N. N. Krasovski, ‘On the Theory of Optimum Control,”’ Prik/. Mat. Mekh., Vol. 23, 
No. 4, pp. 625-639 (1959). 
E. B. Lrg, ‘Mathematical Aspects of the Synthesis of Linear Minimum Response-Time 
Controllers,” JRE Trans. on Automatic Control, Vol. AC-5, No. 4, pp. 283-289 (1960). 
(10) R. E. Katman, ‘Optimal Nonlinear Control of Saturating Systems by Intermittent 
Action,” 1957 IRE WESCON Convention Record, Part 4, pp. 130-135. 
(11) C. A. DESOER AND J. WING, ‘‘An Optimal Strategy for a Saturating Sampled Data Sys- 
tem,’’ JRE Trans. on Automatic Control, Vol. AC-6, pp. 5-15, February, 1961. 
(12) C. A. DESoER AND J. WrnG, “A Minimal Time Discrete System,’’ JRE Trans. on Automatic 
Control, Vol. AC-6, pp. 111-125, May, 1961. 
(13) R. E. Kaman, “On the General Theory of Control Systems,’’ Proc. First International 
Congress of Automatic Control, Moscow, 1960. 
(14) B. FrrEpMAN, “Principles and Techniques of Applied Mathematics,’’ New York, John 
Wiley & Sons, 1956. 
(15) T. M. ApostoL, ‘‘Mathematical Analysis, A Modern Approach to Advanced Calculus,” 
Reading (Mass.), Addison Wesley Publishing Co., 1957. 
(16) E. A. Coppincton anp N. Levinson, ‘Theory of Ordinary Differential Equations,” 
New York, McGraw-Hill Book Co., 1955. 
(17) F. R. GAntmacuer, “The Theory of Matrices,"” New York, Chelsea Publishing Co., 1959. 
(18) D. BLACKWELL AND M. A. Girsuick, ‘Theory of Games and Statistical Decisions,” 
New York, John Wiley & Sons, 1954. 


[J. F. 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 


NATIONAL BUREAU OF STANDARDS NEWS 


A NEW MICROCITE MACHINE 


The National Bureau of Standards has recently developed an electro- 
mechanical searching device that has significant advantages for informa- 
tion retrieval in large collections of documents. Developed by Joshua 
Stern of the Bureau's staff, the machine is based on the ‘“‘microcite”’ 
concept! and employs a special punched-card technique known as 
“peek-a-boo.” The simplicity of the design, together with various 
automatic features, makes possible decentralized search facilities so 
that the inquiring researcher may perform the search himself. 

The improved microcite machine was developed, under the sponsor- 
ship of the Department of Defense, to aid the Bureau’s instrument refer- 
ence service in making information on instruments and control systems 
more readily available to Government and industry. 

The problem of effective retrieval of scientific and technological in- 
formation has increased in importance and in difficulty in recent years. 
Much of the current work in this area is directed towards development 
of complex machines that will search a vast memory automatically and 
furnish the documents sought in a desired form. The high cost of such 
machines requires that their operation be centralized and inhibits direct 
contact of the inquirer with the search in progress. The element of 
feedback, which can contribute much to the effectiveness of a search, is 
difficult to incorporate in such centralized operation. 

The importance of feedback in the search process is illustrated by 
the familiar game of ‘“‘twenty questions.’”” A few moments of considera- 
tion makes it clear that the great effectiveness of twenty questions in 
pinpointing information stems principally from the fact that each 
question is phrased after the answer to the previous question has been 
provided. This same advantage is inherent in the microcite machine. 
Here the researcher is at all times in control of the search process and is 
immediately informed of the nature of the information being discovered. 

A brief description of the peek-a-boo technique and the microcite 
concept should aid in understanding the present development. The 
peek-a-boo technique uses punched cards each of which represents an 
index word. The positions of holes in the card identify documents per- 
taining to the index word. One search techique often referred to as the 
logical product is performed by superimposing the appropriate peak-a- 


1“Microcite, an Aid to More Effective Referencing,’’ NBS Tech. News Bull., Vol. 41, 
p. 141 (1957). 
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boo cards and noting the positions of holes that remain unobscured by 
any of the selected cards. More complex logical operations can also be 
performed. In the ordinary use of the peek-a-boo system, the position 
of a hole is interpreted as a document serial number by means of a co- 
ordinate grid printed upon the card or upon a transparent overlay. The 
searcher is directed by these serial numbers to a file (serial list) of cita- 
tions and abstracts of the documents sought or to a file of the documents 
themselves. 

The microcite concept eliminates this latter step by immediately 
providing the searcher with a microfilm image of the citation and ab- 
stract. Instead of reading serial numbers, the searcher views directly 
at each hole position a description of a document emerging from the 
search. Ina preliminary model, microimages of citations and abstracts 
were reproduced on a film overlay and read at the holes with the aid of 
a microscope. Such a procedure demonstrated the effectiveness of the 
microcite principle but would be relatively inefficient for large collections. 

The new microcite is a tool that could search millions of documents. 
In its particular application at the Bureau, it has been designed for sets 
of 18,000-document peek-a-boo cards with corresponding abstracts for 
each set carried in a 15-in.-square area of a sheet of photographic film— 
the matrix. In operation, the film matrix appropriate to the peek-a-boo 
set being searched is selected and mounted on a drum. The operator 
places selected peek-a-boo cards on an illuminated area of the machine 
and manipulates two control wheels to position a set of cursors on each 
illuminated—that is, unobscured—hole in turn. Turning the control 
wheels causes rotation and translation of the drum carrying the film 
matrix, so that as the cursor is set on each hole, a full-size focused image 
of the citation and abstract corresponding to that hole is projected on a 
screen in front of the searcher. 

As he proceeds, the searcher can obtain a copy of any of the abstracts 
he desires by moving a lever to record the desired image on a roll of 
photosensitive paper. Alternatively, he may note the serial numbers of 
the documents or may cause the serial numbers to be transcribed into 
punched tape or printed on paper tape. The electrically coded serial- 
number output that will be available for this transcription could also 
be fed into a document-storage machine—such as are expected to become 
available in the near future—which would provide copies of complete 
documents in response to serial-number input. 

The immediate availability of the results is a powerful aid to the 
searcher in allowing him to guide the search effectively. He can rapidly 
detect and adjust for such contingencies as misinterpretation of index 
terms by indexer and searcher, unanticipated interference from in- 
formation of whose existence he was unaware, misconceptions about the 
volume of relevant information in the collection, and even modification 
of his interests as he acquires more and more information. 
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POLYMER CRYSTALS WITH FOLDED CHAINS 


The nature of polymer crystals is of considerable interest, since such 
crystals affect the mechanical, optical, and dielectric properties of these 
useful materials. As part of a long-range program involving the dielec- 
tric properties of polymers, scientists at the National Bureau of Stand- 
ards have investigated the fundamental characteristics? of these ma- 
terials from a theoretical point of view. 

In 1958, it was discovered that certain polymers, such as poly- 
ethylene, are deposited from dilute solution in the form of flat and very 
thin crystalline platelets. (See, for example, the work of A. Keller 
and co-workers.*:*) The surprising feature of these plate-like crystals is 
that they are built up of folded polymer chains. The folds themselves 
actually form the large upper and lower surfaces of the platelets. The 
chain axes of the unfolded portions of the polymer molecules are roughly 
perpendicular to the plane of the chain folds. This is a vastly different 
structure from the bundle-like or sheaf-like structure that one might 
have expected the crystals to exhibit. 

A detailed theoretical interpretation of this phenomenon has recently 
been presented by J. I. Lauritzen and J. D. Hoffman of the Bureau’s 
Dielectrics Section. Their explanation involves nucleation and growth 
theory. They have demonstrated that in sufficiently dilute solution, 
the rate of formation of the seeds leading to folded crystals must exceed 
that of seeds leading to crystals of the bundle-like type. 

The growth of each folded nucleus was then considered. It was 
shown that the thickness (often called the step height, /*) will not in- 
crease rapidly as the platelet becomes larger in the x- and y-directions, 
thus preserving the thin platelet character of the crystal. This tendency 
for the step height to maintain itself during the growth process is a 
result of the fact that a crystal of step height /* grows much faster than 
one that is either substantially thicker or thinner. The quantity /* 
generally ranges from roughly 100 to 500 A, and increases with tempera- 
ture of crystallization. 

The theoretical study has recently been extended to the case of crys- 
tallization in bulk polymers.’ Results were obtained which strongly 
suggest that substantially chain-folded crystals (lamellae) exist in bulk 
polymers, especially in the spherulites. Calculations have recently been 
made at the Bureau which give information on the smoothness of the 
chain-folded surfaces of the lamellae and on the details of the melting 
behavior to be expected for chain-folded crystals. 


2 For further information, see ‘Theory of Formation of Polymer Crystals with Folded 
Chains in Dilute Solution,” by John I. Lauritzen, Jr., and John D. Hoffman, NBS J. Research, 
Vol. 64A (Phys. & Chem.), No. 1, p. 73 (1960). 

3A. KELLER AND A. O'Connor, Discussions, Faraday Soc., Vol. 25, p. 114 (1958). 

4A. Kevver, Phil. Mag., Vol. 2, p. 1171 (1957). 

5 “Theory of rate of Lamellar Spherulitic Crystallization in Bulk Polymers,”’ by John D. 
Hoffman and J. I. Lauritzen, Jr.. NBS J. Research, Vol. 65A (Phys. & Chem.), No. 4, p. 297 


(1961). 
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THE FRANKLIN INSTITUTE 


NEW BOOKS IN THE FRANKLIN INSTITUTE LIBRARY 


HISTORY OF SCIENCE AND TECHNOLOGY 


NEWTONIAN SCIENCE, by E. A. Bell. E. Arnold, 1961. 
A SHort History OF TECHNOLOGY FROM THE EARLIEST TIMES TO A.D. 1900, by T. K. Derry. 
Oxford University Press, 1961. 
ROCKET DEVELOPMENT, edited by Esther C. Goddard. Prenctice-Hall, 1961. 
THE NATIONAL WATERCRAFT COLLECTION, by H. I. Chapelle. Smithsonian Institution, 1960. 
DiscouRSE ON BopIEs IN WATER, by G. Galilei, translated by T. Salusbury. University of 
Illinois Press, 1960. 
ASTRONOMY 


INTERNATIONAL SPACE SCIENCE Symposium, edited by H. K. Bijl. Interscience, 1960. 

Space Astropuysics, by W. Liller. McGraw-Hill, 1961. 

TOOLS OF THE ASTRONOMER, by G. R. Miczaika and W. M. Sinton. Harvard University 
Press, 1961. 

SCIENCE IN SPACE, edited by L. V. Berkner. McGraw-Hill, 1961. 


MATHEMATICS 


TRANSFORM CALCULUS FOR ELECTRICAL ENGINEERS, by R. Légros. Prentice-Hall, 1961. 

CONCEPTS FROM TENSOR ANALYSIS AND DIFFERENTIAL GEOMETRY, by T. Y. Thomas. 
Academic, 1961. 

LECTURES ON LINEAR ALGEBRA, by I. Gel’fand. Interscience, 1961. 

ELEMENTS OF MATHEMATICS, by J. H. Banks. 2nd ed. Allyn and Bacon, 1961. 

PARTIAL DIFFERENTIAL EQUATIONS AND CONTINUUM MECHANICS, by R. E. Langer. Univer- 
sity of Wisconsin Press, 1961. 

THEORY OF MARKOV Processes, by E. B. Dynkin. Pergamon, 1960. 

FOURIER TRANSFORMS, by R. Goldberg. Cambridge University Press, 1961. 


PHYSICS 


SMALL PaRTICLE Statistics, by G. Herdan. 2nd ed. Academic, 1960. 

TURBULENCE; CLASSIC PAPERS ON STATISTICAL THEORY, edited by S. K. Friedlander. I nter- 
science, 1961. 

DIFFRACTION, STRUCTURE DES IMAGES; INFLUENCE DE LA COHERENCE DE LA LUMIERE, by 
A. Maréchal. Editions de la Revue d'optique théorique et instrumentale, 1960. 

HYDROMETRY; THEORY AND PRACTICE OF HYDRAULIC MEASUREMENTS, by A. T. Troskolafiski. 
Pergamon, 1960. 

SYMPOSIUM ON RADIATION EFFECTS AND RADIATION DosIMETRY. American Socity for Testing 
Materials, 1960. 

THE FUNDAMENTAL Atomic Constants, by J. H. Sanders. Oxford University Press, 1961. 

ELECTROMAGNETIC STRUCTURE OF NuCLEONS, by S. D. Drell. Oxford University Press, 1961. 

HEAT TRANSFER, by B. Gebhart. McGraw-Hill, 1961. 

ELEMENTS OF NUCLEAR ENGINEERING, by G. Murphy. Wiley, 1961. 

ATOMIC ENERGY WASTE, by E. Glueckauf. Interscience, 1961. 

Puysics AND ARCHAEOLOGY, by M. J. Aitken. Interscience, 1961. 

PHYSICS FOR ENGINEERS AND SCIENTISTS, by R. G. Fowler. Allyn and Bacon, 1961. 

THE PRINCIPLES OF NUCLEAR MAGNETISM, by A. Abragam. Clarendon Press, 1961. 
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CHEMISTRY AND CHEMICAL ENGINEERING 


REINFORCED PLASTICS FOR ROCKETS AND AIRCRAFT. Symposium. American Society for 
Testing Materials, 1961. 

INORGANIC CHEMISTRY, A GUIDE TO ADVANCED Stupy, by R. B. Heslop. Elsevier, 1960. 

INTERNATIONAL SYMPOSIUM ON MACROMOLECULAR CHEMISTRY.  Izdatel’stvo Akademii Nauk 
SSSR, 1960. 

PoLYMERIC MATERIALS, by C. C. Winding. McGraw-Hill, 1961. 

Cuemistry, by M. J. Sienko. McGraw-Hill, 1961. 

INTRODUCTION TO QUALITATIVE ANALYsIS, by D. C. Layde. Allyn and Bacon, 1961. 

REAGENT CHEMICALS; SPECIFICATIONS, 1960. American Chemical Society, 1961. 

NAME INDEX OF ORGANIC REACTIONS, by T. S. Wheeler. Interscience, 1960. 

THE SCIENCE OF ADHESIVE JOINTS, by J. J. Bickerman. Academic, 1961. 

SYMPOSIUM ON ADHESION AND ADHESIVES. American Society for Testing Materials, 1961. 

KOLLOIDCHEMISCHES TASCHENBUCH, by A. Kuhn. Akademische Verlagsgesellschaft, 1960. 

INFRARED ABSORPTION OF INORGANIC SUBSTANCES, by K. Lawson. Reinhold, 1961. 


MECHANICAL ENGINEERING 


CREATIVE ENGINEERING ANALYsIS, by F. L. Ryder. Prentice-Hall, 1961. 

THEORIES OF ENGINEERING EXPERIMENTATION, by H. Schenck. McGraw-Hill, 1961. 

THe CONSULTING ENGINEER, by C. M. Stanley. Wiley, 1961. 

ENGINEERING DATA FOR PROpucT DersIGN, by D. C. Greenwood. McGraw-Hill, 1961. 

50 YEARS PROGRESS IN MANAGEMENT, 1910-1960. American Society of Mechanical Engineers, 
1960. 

SMALL PLANT MANAGEMENT, by W. A. MacCrehan, Jr. 2nd rev. ed. McGraw-Hill, 1960. 

MECHANICAL MEASUREMENTS, by T. G. Beckwith. Addison-Wesley, 1961. 

HANDBOOK OF MECHANICAL WEAR, by C. Lipson. University of Michigan Press, 1961. 

PLastic by M. W. Riley. 2nd ed. Reinhold, 1961. 

FRICTION AND WEAR IN MACHINERY. American Society of Mechanical Engineers, 1960. 

Tue ROLe or Viscosity IN LUBRICATION, edited by O. C. Bridgeman. American Society of 
Mechanical Engineers, 1960. 

DICTIONARY OF MECHANICAL ENGINEERING, by A. Del Vecchio. Philosophical Library, 1961. 

THe DESIGN OF SMALL DiREcT-CURRENT Morors, by A. F. Puchstein. Wiley, 1961. 

MECHANIZATION OF Motion, by L. Harrisburger. Wiley, 1961. 


ELECTRICAL ENGINEERING 


ELECTRICAL ENGINEERING, by S. B. Hammond. McGraw-Hill, 1961. 

AN INTRODUCTION TO THE PRINCIPLES OF COMMUNICATION THEORY, by J. C. Hancock. 
McGraw-Hill, 1961. 

COMPUTER ENGINEERING, edited by S. A. Lebedev. 


Pergamon, 1960. 


METALLURGY 


SymposiuM ON Low-TEMPERATURE PROPERTIES OF HIGH-STRENGTH AIRCRAFT AND MISSILE 
MATERIALS. American Society for Testing Materials, 1961. 

FATIGUE TESTING AND ANALYSIS OF RESULTS, by W. Weibull. Pergamon, 1961. 

INTERNATIONAL SYMPOSIUM ON THE PHYSICAL CHEMISTRY OF PROCESS METALLURGY, edited 
by G. R. St. Pierre. Interscience, 1961. 

ENGINEERED CAstTINGs, by G. J. Cook. McGraw-Hill, 1961. 

EXTRACTIVE METALLURGY OF CopPrpER, NICKEL, AND COBALT, edited by P. Queneau. Inter- 


science, 1961. 
CosaLt MonoGrapn, prepared in collaboration with the staff of Battelle Memorial Institute. 


Brussels, 1960. 
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BOOK REVIEWS 


AERO- THERMODYNAMICS AND FLOW IN TURBO- 
MACHINES, by Michael H. Vavra. 609 
pages, diagrams, 5 X 9 in. New York, 
John Wiley & Sons, Inc., 1960. Price, 
$14.50. 


This book presents a vector and dyadics 
analysis of fluid flow through turbomachinery. 
It includes applications to axial flow compres- 
sors and turbines and two-dimensional poten- 
tial flow through blading cascades. <A sim- 
plified three-dimensional analysis of axial 
turbomachines is included. Other features 
include an extensive appendix on vector 
analysis, including an extension to operational 
analysis involving dyadics. 

The book is divided into three major parts. 
Part I pertains to ‘Fundamental Flow Rela- 
tions’; Part II extends to relative and rota- 
tional motion, and also discusses isentropic 
and potential flows; Part III covers applica- 
tions to turbomachines already mentioned. 

In Part I, the author formulates the thermo- 
dynamic properties in the equations of motion, 
where the latter are expressed in terms of 
gradients of velocity and pressure. Inserting 
the thermodynamic property 


1 
T-As = Vh —-Vb 
for the pressure drops in 
1 
—+VVV = --VWp-F, 
ot e 


the author arrives at 


OV = —vh4 TVs — F, 
at 

where the fluid flow can be expressed in terms 
of enthalpy and entropy. Since the increase 
of entropy, As, is the combined effect of 
both external and internal heat transfers, 
while the irreversibility due to friction is it- 
self TAsr = fdF, approximately, the ad- 
vantage of this method appears only by ap- 
proximating entropy increases, but retaining 
explicitly the friction forces f. Actually, as 
with steady flow, with adiabatic irreversible 
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V-VV = — Wh, due to the above 
relation. The irreversibility can likewise be 
accounted for by comparing actual with 


systems, 


isentropic enthalpy changes. 

The arguments relating to energy relations 
in the author's discussion of jet propulsion 
seem highly artificial, if not incorrect. In 
estimating the jet reaction, the velocities 
V; and V2 can be interpreted either as rela- 
tive or velocities, for the 
“change” in momentum, we are only con- 
cerned with the difference in velocities. How- 
ever, in the statement of energy relations, it 
is very important to distinguish between the 


absolute since 


two. 


The author (p. 75) states H; = hy a4 
g 


7,2 
and Hz = hz + — as the total energies at 


stations 1 and 2, where V; and Vs are in- 
terpreted as relative velocities. With the 
unnecessary hypothesis to prove a funda- 
mental relation, the author assumes an “‘ideal”’ 
engine with enthalpies h; = hz and T; = T2 


and concludes — H, = 


minimum possible energy for acceleration of 
a particle from V; to V2. It is important to 
dM, 

point out that f 

V; and V2 interpreted as relative velocities, 
is always a measure of the ‘mechanical’ 
energy imparted to the flow stream, irrespec- 
tive of whether the thermal cycle is reversible 
or irreversible. AE should be regarded as the 
available part of the heat input, AQ, for 
mechanical work. Hence any assumption 
as to the nature of the thermal cycle is 
unnecessary. 

Actually, the total energy for mechanical 
work must include the heat energy injected 
into the flow stream through combustion of 
the fuel supply at a rate AQ per unit time. 
The heat energy retained in the flow stream 
at exit is measured by the exit enthalpy. 
With the convected enthalpies 4; and he at 
terminal sections, a flow rate G(lb/sec), w as 
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relative, and v = the velocity of plane, then 
AE is 


G 
AE = AQ + — hz) = (we? — 
g 


which gives the relation of entrance and exit 
enthalpies, with heat source AQ. Neglecting 
the small mass flow of injected fuel, it can be 


G 
shown that ~ (ws? — w,*), that is, the square 


of the difference of the relative velocities, ac- 
counts for the heat energy transferred to 
mechanical work, part of which is utilized 
in propelling the plane and the other is lost 


in the slip stream. Even with the author’s 


G 
ideal engine, if AE = (w.? — w,?) = AQ, 


with h; = he, then 7; = T2 must be nil. The 
author's ideal engine appears inconsistent with 
the Carnot cycle, where only a part of AQ 
is available for mechanical work, if finite 
exit temperatures are assumed where hz > hy. 
A possible confusion arises in not clearly 
separating the heat energy transferred to 
mechanical work, and the heat energy retained 
in the slip stream and finally transferred to 
the atmosphere at its ambient temperature. 
This phase is concerned with the thermal 


efficiency. The energy transformed to me- 
chanical work, again separates to propulsion 
work and to mechanical energy lost in the 
slip stream. This latter is related to the 
propulsion efficiency which the author is 


concerned with. The propulsion efficiency 


2u 
€p = — w)v/(w2? — w,?) = ————,, 
(we + w1) 
where w,; = v = velocity of plane. On the 


other hand, the thermal efficiency is 


— wi?) 


AQ 
and the over-all efficiency is 
G (we — 

AQ 
Moreover, whether the energy transferred is 
isentropic or irreversible, the previous equa- 
tion and the above relations hold. Therefore, 
the author’s hypothetical ideal cycle appears 


unwarranted for evaluating the propulsion 
efficiencies. 
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In Part II, the author discusses relative 
motions of the flow through rotors. The 
distinction between steady flow in stators 
and unsteady flow in rotors with respect to 
stator, is clearly brought out. With steady 
relative flow in rotor, the significance of 
rotational motion, the Coriolis and centrif- 
ugal component and circulation flow, are dis- 
cussed. The significance of the Coriolis com- 
ponent and the rate of change of angular 
momentum should be emphasized. The char- 
acteristics of rotational flow, with vortex 
filaments are mathematically discussed. In 
this section isentropic flow of perfect gases 
is discussed, followed by potential flow 
analysis. 

The remainder of the book includes applica- 
tions to axial flow compressors and the analy- 
sis of cascade flow. 

Special mention should be given to the 
author's very complete review of vector analy- 
sis, and an unique introduction to operations 
involving dyadics. 

This book is both an interesting and valu- 
able discussion of fluid flow analysis in turbo- 
machinery. It extends a type of mathematics 
that seems indispensable to advancing the 
analysis and research studies of fluid flow 
problems. 

RuPrEN EKSERGIAN 
The Franklin Institute Laboratories 


ADVANCES IN THE ASTRONAUTICAL SCIENCES, 

VoL. 6, edited by Horace Jacobs and Eric 
898 pages, diagrams, illustra- 
New York, The Mac- 
Price, $25.00. 


Burgess. 
tions, 63 in. 
millan Co., 1961. 


The annual meetings of astronautical so- 
cieties permit the presentation of papers which 
mirror both current events and thinking in the 
space-travel panorama. These meetings allow 
leaders in allied fields to explore facets of 
astronautics of common interest. Thus at 
these sessions can be heard discussions of the 
most important topics of the day. 

In January 1960 the Sixth Annual Meeting 
of the American Astronautical Society was 
held in New York. The three-day session 
was the setting for 56 technical papers, three 
special lectures, and a space panel. This 
volume presents these papers. 

There are nine sections to the book which 


‘ 
: 
3 
: 


236 


range from communications to space physics. 
Some of the sections, such as Guidance and 
Control and Space Flight Mechanics, are 
highly developed technical areas which need 
a substantial background in the subjects to 
understand them. Other sections are of a 
more popular nature. 

One of the papers in the space-flight 
mechanics section by E. T. Benedikt 
cerns “Disintegration Barriers to Extremely 
High-Speed Space Travel.’’ The author’s 
thesis is that when we travel at high enough 


con- 


velocities, the space ship will disintegrate 
due to the impacts of particles we meet in 
space. In the reviewer's book “Target For 
Tomorrow” this writer showed that indeed 
this was a possibility and the only way to 
circumvent the hazard was to fabricate the 
interstellar vehicle from an asteroid leaving 
a shell several hundred feet thick. Thus 
while the material would wear away there 
would still be enough material there to main- 
tain the integrity of the space vehicle. 
Dandridge M. Cole contributed a paper on 
“The Feasibility of Propelling Vehicles by 
Contained Nuclear Explosions."” This re- 
viewer has always admired the boldness of 
this concept and here it is detailed for the 
reader. Mr. Cole in his paper is concerned 
with all-up weights of 175,000 tons 
than twice the weight of the largest ocean 


more 


liner. 
siderations of Expected Radiation Belts of 
Planets Mars and Venus,” is another example 
of the analytical work being done by a rank- 
Dr. 
Singer may be remembered as the scientist 
who in 1958 predicted and explained the radia- 
tion belts of the earth before their discovery. 


Fred Singer’s paper, “Some Con- 


ing scientist in the space-travel field. 


Here is another example of his ingenious de- 
ductions. Now it remains for us to determine 
experimentally whether he is right in his 
analysis. 

These are but a few of the papers in this 
outstanding volume which caught the interest 
of this reviewer. Everyone in the space- 
travel field will find papers of extreme interest 
to him. 

Advances in the Astronautical Sciences is 
another reference book necessary to permit 
one to stay abreast of this exploding science. 

I. M. Levitt 
The Fels Planetarium 
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ELECTROMAGNETIC FIELDS AND WaAvVEs, by 
R. V. Langmuir. 227 pages, diagrams, 
6 X9in. New York, McGraw-Hill Book 
Co., Inc., 1961. Price, $9.75. 


Professor Langmuir states in this new book 
that “In such a rapidly moving field as modern 
electrical engineering one can scarcely do 
more than train the student in the 
mathematical methods and basic physics of 
the field . . . to meet and solve problems 
that now.”” He 
expects the students to come quite well pre- 


much 


are not even dreamed of 
pared in vector analysis, partial differential 
Fourier Thus, his 
new rather little 
time is spent tutorially in building mathe- 
matical background or in discussing current 


series, etc. 


short, 


equations, 


text is because 


practices. 

What is covered is treated in a precise, ele- 
gant, and sophisticated way. It is safe to 
say that the quality of the book rests partly 
with the author’s good judgment of what to 
leave out. Of course there is no new material; 
the entire text, having only 200 pages, cer- 
tainly had to be prepared with strict dis- 
crimination. For example, the treatment of 
waveguides does not contain data on the nu- 
merous guide configurations and irises. How- 
ever, the exposition is sufficiently powerful to 
arm the student with the means for analyzing 
a great many of the microwave engineering 
problems which may be encountered. 

Advanced subjects such as plasma physics, 
about which engineers of the future will need 
to know much, are not covered at all. On 
the other hand, entire chapters are devoted 
to the fundamental Laplace and Maxwell 
equations. The mathematical tech- 
niques for dealing with dielectric phenom- 
ena and magnetic materials are similarly 


basic 


treated. 

In summary, it is apparent that, through 
the courses he has taught at the California 
Institute of Technology, Professor Langmuir 
has distilled certain basic elements from such 
sources as Smythe, Stratton, Slater, Frank, 
et al. In teaching this highly mathematical 
subject the author makes a point of empha- 
sizing the tangible reality of terms by the use 
of some delightful, homespun examples such 
newton about the gravitational 
force on a small apple.’’ The result is an 
excellent, compact and concentrated version 


as “3 
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of the classic texts for effective use by the 
senior electrical engineering student. 

C. W. HAarRGENS 

The Franklin Institute Laboratories 


MEN OF SpacE—VOLUME II, by Shirley 
Thomas. Chilton Press, 1961. Price, 
$3.95. 


For Volume II of her series on the men who 
are shaping our destinies in space, Shirley 
Thomas has chosen ten more distinguished 
scientists. 

While Volume I profiled the obvious first 
choices, with the names of the men who 
figured most prominently in the news, this 
volume is by no means a collection of bio- 
graphical sketches of the second team. We 
must face up to the situation that in this 
country the news media can promote and 
glamorize an individual who may say and do 
no more or less than a close associate or col- 
league. We never hear of that colleague, yet 
he is an important contributor to space 
science. Thus this volume is essentially one 
detailing the contributions of the lesser known 
but still brilliant stars in the space firmament. 

Some of the men have obvious connections 
with space as in the case of Dr. Dornberger, 
Dr. Whipple, Captain Truax, Dr. Lovelace, 
etc. Others, Drs. Ramo, Teller, 
though dynamic figures, are involved only 
indirectly in the space-flight picture. 

The author has a unique ability of portray- 
ing the subject in vignettes which spotlight 
significant events. Thus the reader is per- 
mitted to fill in the obvious details. By 
writing in such fashion, the 225-page book 
represents a much more complete 
She has the happy faculty of por- 


such as 


really 
volume. 


traying the scientist not only as a worker in 
his discipline but as a well-rounded human 
being who has a family and is subject to all 
the problems and situations which arise from 
this 


state. Biographies make fascinating 
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reading and the author has the enviable 
faculty of making her subjects come alive. 

The significance of these Men in Space 
profiles will be appreciated later when their 
pioneering work becomes history. Our chil- 
dren will be better off for knowing just how 
this adventure started and who gave it the 
greatest push and how they did it. 

Miss Thomas has again produced a live, 
stimulating and highly readable book. I look 
forward to the next in the series. 

I. M. Levitr 
The Fels Planetarium 


Astropaysics, edited by William 

272 pages, diagrams, illustrations, 
6X9 in. New York, McGraw-Hill Book 
Co., Inc., 1961. Price, $10.00. 


During the 1959-1960 year the Astronomy 
Department of the University of Michigan 
presented a series of lectures on those aspects 
of astronomy and astrophysics which could be 
investigated only from outside the atmosphere 
of the earth. 

This volume is a compilation of the lectures 
given by ranking authorities at the invitation 
of the University. Their contributions are 
in most cases of original material and find 
their way into print for the first time. 

The book is divided into two sections with 
the first section involving the experiments 
performed and to be performed in space, 
and the second section is concerned with the 
theory. The second section is an excellent 
“wrap-up” of our knowledge of the solar 
wind, corona, the electric universe and the 
theory behind some of the instruments which 
may be used to further our knowledge of 
space. 

The second section makes it a worth-while 
addition to the library of the working space 


scientists. 


SPACE 
Liller. 


I. M. Levitt 
The Fels Planetarium 


ae 
‘Gua 
x 
ak 
: 
: 


BOOK NOTES 


PROGRESS IN BIOPHYSICS AND BIOPHYSICAL 
CHEMISTRY, VOL. 11, edited by J. A. V. 
Butler, B. Katz and R. E. Zirkle. 255 
pages, diagrams, plates, 6} 10 in. Lon- 
don, Pergamon Press, 1961. Price, $12.50. 


The overlapping of the physical and the 
biophysical sciences is emphasized in this 
series. Volume 11 is comprised of six papers, 
with an added section of short discussions 
covering various phases of cytoplasmic par- 
ticles and their role in protein synthesis. The 
six major papers are: ‘The Natural Radio- 
activity of the Human Body”; Engi- 
neering Approach to the Problem of Neural 
Organization” ; “Effects of X-Rays on Nucleic 
Acid Biosynthesis and on the Activity of 
Nucleases in Mammalian Cells’’; ‘“The Effect 
of Ionizing Radiations and Tumour-Chemo- 
therapeutic Agents on the Bone Marrow’’; 
“Some Factors Influencing the Dispersion of 
Indicator Substances in the Mammalian 
Circulation’; and ““The Buoyancy of Fish 
and Cephalopods.”’ Excellent bibliographies, 
including 1960 references, accompany not 
only the six papers, but also the seventh sec- 
tion on cytoplasmic particles. 


” 


ENERGY PRINCIPLES IN APPLIED STATICS, by 


T. M. Charlton. 106 pages, diagrams, 
6X9 in. London, Blackie & Son Ltd., 
1959. Price, 25s. 


The author has written this text in an effort 
to clarify the difference between the strain- 
energy and potential-energy methods appli- 
cable to the analysis of statical, plane pin- 
jointed systems. The author clearly demon- 
strates the principle of virtual work, and 
points out that Castigliano’s “principle of 
least work”’ (which is a misnomer) is only one 
of two possible approaches to the analysis of 
linear statically-indeterminate structures, the 
other being the method of potential energy. 
The author devotes two (of a total of seven) 
chapters to the complementary-energy method 
of analysis. Our readers may recall Charlton's 
paper in this JouRNAL (Vol. 250, p. 543) on 
the analysis of redundant systems through 
the conception of conservation of energy. 
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DISPERSION RELATIONS & ELEMENTARY PAr- 
TICLES, edited by Cécile De Witt and R. 
Omnes. 671 pages, diagrams, 6} X 9} in. 
Paris, Hermann & Cie., 1961. Distributed 
in U.S. by John Wiley & Sons, Inc. Price, 
$20.00. 


This collection of eight treatises on various 
aspects of dispersion relations and elementary 
particles represents material presented at the 
1960 summer session of the School of Theo- 
retical Physics at the University of Grenoble. 
Two of the papers are in French, the others in 
English. 

M. L. Goldberger’s “Introduction to the 
Theory and Applications of Dispersion Rela- 
tions”’ covers the historical background, two- 
dimensional representation, non-relativistic 
and relativistic theory, and applications to 
pion-nucleon scattering, forward scattering 
and non-forward scattering. 

A. S. Wightman’s two papers (one in 
English and one in French) deal with invari- 
ance in relativistic quantum mechanics (Lor- 
entz and Poincaré groups, physically and not 
physically realizable states, inversions, in- 
finite operators) and with analytic functions 
of several complex variables (holomorphy 
domains, continuity theorem, and others). 

R. Omnes’ contribution (in French) covers 
general properties of the S matrix reduction 
formulas, Dyson's formula, and holomorphy 
domains, with mathematical proofs and 
examples. 

A. O. G. Kallen discusses the properties of 
vacuum expectation values of field operators, 
covering field and particle concepts, reduction 
formulas, 2-, 3-, 4- and n-point functions. 

G. F. Chew's contribution, ‘‘Double Dis- 
persion Relations and Unitarity as the Basis 
for a Dynamical Theory of Strong Interac- 
tions,”’ deals, for example, with the Lorentz- 
invariant amplitude, the Mandelstam repre- 
sentation, charge and spin, partial-wave am- 
plitudes, pion-pion scattering, and the nucleon 
electromagnetic structure. 

S. B. Treiman discusses weak interactions, 
including Leptonic and non-Leptonic reac- 


tions, and universal models. Included are 
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beta decay, u-meson decay, decay of neutral 
K particles, and others. 

The final contribution, by Y. Yamaguchi, 
deals with interactions of strange 
particles, covering isospin, spin-parity, K-N 
collisions, the Stark effect, and the application 
of dispersion theory to polology. 


strong 


PRINCIPLES OF INTERNAL NAVIGATION, by 
C. J. Savant, Jr., R. C. Howard, C. B. 
Solloway and C. A. Savant. 254 pages, 
illustrations, 6X9 in. New York, 
McGraw-Hill Book Co., Inc., 1961. Price, 
$9.75. 

This first book devoted entirely to this 
rapidly expanding field covers the basic con- 
cepts of automatic navigation, the design and 
testing of autonavigation components, and 
the necessary mathematical analysis tech- 
niques. Recent advances are included, such 
as ballistic guidance and aircraft guidance. 
The text, which includes most of the unclassi- 
fied material available, was originally prepared 
as a series of lectures for graduate students, 
but it is also a valuable reference for those 
working in the field. Among the 47 refer- 
ences (dated 1954-1959) O'Donnell 
paper appearing in this JOURNAL for 1958. 


is the 


ZAHLENWERTE UND FUNKTIONEN AUS PHYSIK, 
CHEMIE UND ASTRONOMIE, GEOPHYSIK UND 
Tecunik. II Banp: EIGENSCHAFTEN DER 
MATERIE IN IHREN AGGREGATZUSTANDEN, 
7 Tet: ELEKTRISCHE EIGENSCHAFTEN II 
(ELEKTROCHEMISCHE SYSTEME), by R. Ap- 
pel, K. Cruse, P. Drossbach, H. Falken- 
hagen, G. G. Grau, G. Kelbg, E. Schmutzer 
and H. Strehlow, edited by Karl-Heinz 
Hellwege, Anne Marie Hellwege, Klaus 
Schaffer and Ellen Lax. Sixth edition, 
959 pages, diagrams, 74 X 10} in. Berlin, 
Springer-Verlag, 1960. Price, DM 478. 


This second section of Part 7, which has 
the general title of Electrical Properties, 
covers electrochemical systems in the careful 
and detailed manner with which all the vol- 
umes in this important series of tables and 
functions have been prepared. This sixth 
edition is divided into three main parts— 
Electrical Conductivity in Electrochemical 
Systems, Electromotive Forces, and Equi- 
librium in Electrochemical Systems. The 
first part covers tables of molten salts, pure 
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liquids, aqueous and nonaqueous solutions 
and electrophoretic and electrokinetic poten- 
The second part contains tables of 
The 
third part covers dissociation constants, acids, 
bases and indicators, and pH numbers with 
buffer solutions. The table of contents is 
clearly arranged so that a particular table may 
be quickly located. Literature sections follow 
all of the tables. The book is well printed and 
well bound, to withstand the hard usage de- 
manded of a reference work of this sort. 


tials. 
reversible and irreversible potentials. 


THE PLASMA StaTE, by E. J. Hellund. 190 
pages, diagrams, 6 X 9 in. New York, 
Reinhold Publishing Corp., 1961. Price, 
$6.50. 


This book is an initial effort in presenting 
the fundamentals and significance of highly 
ionized gases, fully describing what is cur- 
rently known of the plasma state in all its 
aspects. 

The author defines plasma and explains the 
complexities of energy-exchange processes, 
after giving a brief historical outline of high- 
energy sources. Plasma morphology, energy 
exchange and hot gases are covered in separate 
chapters. A discussion of problems yet to be 
solved in regard to plasma is included in a 
chapter on the various kinds of plasmas, and 
there is an outline of devices, and information 
on the artificial and natural sources of high- 
energy particles. Industry research 
workers will find a wealth of valuable informa- 
tion in this important new book. 


and 


TopoLocy, by John G. Hocking and Gail S. 
Young. 363 pages, diagrams, 6 X 9 in. 
Reading (Mass.), Addison-Wesley Pub- 
lishing Co., Inc., 1961. Price, $8.75. 


Offering a unified and self-contained treat- 
ment of basic topology, this book is intended 
for a first course in the subject of twosemesters 
in length, at the graduate level. The authors 
had a double goal in writing this text—to 
present the fundamentals of topology needed 
by every professional mathematician and also 
to erect a broad framework upon which the 
student of topology may build. As a result 
of this, almost every topic of interest in 
topology is mentioned. 

The first half of the book is primarily 
set-theoretic and the second half is algebraic, 
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being particularly noteworthy for its detailed 
treatment of the elementary homotoph and 
homology theory of simplicial complexes. 
Some of the topics covered are the theory of 
continua, relative homology and cohomology, 
Peano continua, inverse limit sequence and 
the Eilenberg-Steelrod axioms for simplicial 
theory. 


MOopeERN DIELECTRIC MATERIALS, edited by 
J. B. Birks. 
6X9 in. 
Inc., 1961. 


246 pages, diagrams, plates, 
New York, 
Price, $7.50. 


Academic Press, 


This book should be of value to engineers, 
physicists and chemists who are interested in 
dielectric materials, whether in their develop- 
ment, manufacture or practical applications. 
Providing an introduction to the chief dielec- 
tric materials, there are chapters on paper, 
hydrocarbon insulating oils, chlorinated hy- 
drocarbons, natural and synthetic rubbers, 
synthetic high polymers, silicones, ceramics, 
glass, mica and micanite, and fibres and 
textiles. 

Since there are so many widely different 
materials involved, a team of specialists, 
each of whom is actively engaged in research 
and development in the field of which he 
writes, has been assembled under the editor- 
ship of Dr. J. B. Birks of the University of 
Manchester. 


RADIOACTIVE WastTEs, edited by John C. 
Collins. 231 pages, diagrams, 6 X 9 in. 
New York, John Wiley & Sons, Inc., 1961. 
Price, $8.00. 


One of the most important contributions of 
recent years to the field of wastes engineering, 
this book covers the implications of radio- 
activity in regard to water supply, waste 
water disposal and disposal of radioactive 
wastes and radioactive gases. Of a funda- 
mental nature, the first half of the book ex- 
plains the problems which have to be under- 
stood when dealing with radioactive material. 
The second half deals specifically with the 
treatment and disposal of liquid, solid and 
gaseous radioactive wastes. This latter part 
also information on 
results of laboratory work and on the experi- 


includes valuable the 


ence gained in operating treatment plants in 
different parts of the world. 
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ELECTRONIC MEASURING INSTRUMENTS, by 
Harold E. Soisson. 338 pages, diagrams, 
illustrations, 6 X 9 in. New York, Mc- 
graw-Hill Book Co., Inc., 1961. Price, $7.50. 


This text will prove an excellent reference 
source for engineers for basic information on 
electronic measuring instruments. pro- 
vides an introduction to electronic measure- 
ment and control as used in the laboratory 
and industry, and offers a general background 
for advanced work. 

The first six chapters acquaint the student 
with the theory and function of the component 
unit, the fundamental 


needed to establish a common background 


covering concepts 
for the laws, terms, basic concepts, symbols 
and nomenclature required for an under- 
standing of the functions and use of electronic 
instruments. The last eight chapters show 
the components’ function in an instrument 
system used to measure or control a labora- 
tory or industrial variable. Some topics of 
special interest which are covered in the book 
are semiconductor diodes and _ transistors, 
radiation counting equipment and analog and 


decimal converter systems. 


PROCEEDINGS OF THE SIXTH CONFERENCE ON 
MAGNETISM AND MAGNETIC MATERIALS, 
sponsored by the American Institute of 
Electrical the American 
Institute of Physics. 395 pages, illustra- 
tions, 8 X 104 in. New York, McGraw- 
Hill Book Co., Inc., 1961. Price, $10.00. 


Engineers and 


Issued as a supplement to the March 1961 
Journal of Applied Physics, this volume con- 
tains the 149 papers presented at the Novem- 
ber 1960 Conference on Magnetism and Mag- 
netic Materials held in New York City. 
Fifteen topics (for example, Ordered Spin 
Systems, Metallic Films, Ferromagnetic Reso- 
nance, Anisotropy, Microwave Devices, etc.) 
are covered, with from five to sixteen papers 
in each topic. Publication of the proceedings 
within a very few months of the conference 
will enhance its value to those in the field. 


TRANSISTOR Locic Circuits, by Richard B. 
Hurley. 363 pages, diagrams, 6 X 9 in. 
New York, John Wiley & Sons, Inc., 1961. 
Price, $10.00. 


The author covers both logical mathematics 
(including logical routines and blocks) and 


ats 
é 
: 
Bee 
i 
: 


Sept., 1961.] 


the transistor circuits that implement them. 
Although developed as a text for seniors and 
graduate students, the book can be profitably 
used by practicing engineers. The book is 
roughly divided equally between switching 
theory and switching circuits. The end-of- 
the-chapter problems are designed to illustrate 
the text material rather than to challenge the 
student to creative thinking. Sample chap- 
ter headings are: Binary Arithmetic; Diode 
Switches and Logic Circuits; Minimization 
(both single function and multiple function) ; 
Sequential Logic—Counters and Registers; 
and Monostable and Astable Circuits. 


HANDBOOK OF FLUID DyNamics, edited by 
Victor L. Streeter. 1240 pages, illustra- 
tions, 6 X 9 in. New York, McGraw-Hill 
Book Co., Inc., 1961. Price, $24.00. 


This is an engineering and scientific level 
handbook, prepared by a panel of more than 
30 well-known engineers and scientists who 
represent a cross section of industry, govern- 
ment laboratories and universities. 

Fluid-flow principles, theory, methods, and 
allied data are broadly covered. Although 
the coverage is highly practical with specific 
applications to the fields of hydraulic power, 
propulsion, aerodynamics, petroleum produc- 
tion, and chemical reaction, theory is included 
whenever necessary for a clearer understand- 
ing of the underlying principles. Divided 
into two parts, the first deals with funda- 
mental concepts and the second is devoted 
to applied fields. 

Many detailed illustrations are used as 
aids for the reader and to provide current 
values of important fluid dynamics data. 


PLAastic AND FRACTURE IN SoLips, by 
Tracy Y. Thomas. 262 pages, diagrams, 
6X9 in. New York, Academic Press, 
Inc., 1961. Price, $8.50. 


Volume 2 of the Mathematics in Science 
and Engineering series, this book will be of 
special use to engineers and students of the 
theory of plasticity in the solution of some 
of their problems. Based largely on material 


from articles published by the author in 
various mathematical journals during the last 
five years, this book is primarily a study of the 
propagation, growth, and decay of discon- 
tinuities in solids. 


The theory of charac- 
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teristic surfaces and the spread of plastic 
into elastic regions is also included. 

The author’s objective discussion of the 
topics meets a proper mathematical standard. 


MEASURE, LEBESQUE INTEGRALS, AND HIL- 
BERT SPACE, by A. N. Kolmogorov and 
S. V. Fomin, translated by Natascha A. 
Brunswick and Alan Jeffrey. 147 pages, 
6X9 in. New York, Academic Press, 
Inc., 1961. Price, $4.00. 


Originally published in Moscow in 1960, 
this is the second volume of the two-volume 
work Elements of the Theory of Functions and 
Functional Analysis. 

The theory of measure and of the Lebesque 
integral, as presented in this text, is based 
on lectures given by A. N. Kolmogorov in the 
Mechanics- Mathematics Faculty of the Mos- 
cow State University. The final preparation 
of this book was carried out by S. V. Fomin. 
The book is divided into five chapters: 
Measure Theory, Measurable Functions, The 
Lebesque Integral, Functions Which Are 
Square Integrable, and The Abstract Hilbert 
Space: Integral Equations with a Symmet- 
ric Kernel. 


RapDIATION BiopHysics, by Howard L. An- 
328 pages, diagrams, illustrations, 
Englewood Cliffs (N. J.), Pren- 
Price, $11.35 (trade) ; 


drews. 

6 X 9 in. 

tice-Hall, Inc., 1961. 

$8.50 (text). 

This new work presents all significant infor- 
mation and pertinent biophysical techniques 
essential to an understanding of radiation 
biology, but never before brought together 
in a textbook. The author, Chief of the 
Radiation Physics Section of the National 
Cancer Institute, provides an important land- 
mark in the study of a subject which has an 
impact on all facets of human existence. He 
presents the workable knowledge and mathe- 
matical background which are needed for 
treating mortality data, analyzing experi- 
mental results for target size and shape, and 
hitness number. 

Including methods of dosage calculation 
not previously available in a textbook, work- 
ing graphs, tables of constants and other data 
required in radiation computation, this book 
will be welcomed by those employed in the 
field of medical science, radiology, biology, 


= 
= 
: 
; 
ae 
: 
Bia 
: 
2 
Re 


242 


radiation biophysics, as well as all those 
concerned with the hazards and the rich 


promise of nuclear power. 


TECHNIQUES OF HIGH ENERGY Puysics, Vol. 
5, edited by David M. Ritson. No. 5 in 
the Interscience Monographs and Texts in 
Physics and Astronomy series. 540 pages, 
diagrams, illustrations, 6 X9 in. New 
York, Interscience Publishers, Inc., 1961. 
Price $16.75. 

This book was written, the editor states in 
his preface, to help bridge the gap ‘“‘between 
the young physicist starting his research 
career and the professional who has grown 
within a given field.””. Among the eleven con- 
tributors to this work are men from IBM 
Corp., M.1I.T., AVCO Research Laboratory, 
and Harvard, Princeton, Brown and North- 
eastern Universities. 

The book is broken 
general headings, which 
Emulsions, Scintillation 
Counters, Beam Optics and General Proper- 
ties of Particles and Radiation. 


eleven 
Nuclear 


Cerenkoy 


down into 
include 


and 


FREE RapicaLs IN BIOLOGICAL SysTEMS, 
edited by M. S. Blois, Jr., H. W. Brown, 
R. M. Lemmon, R. O. Lindblom and M. 
Weissbluth. 387 pages, diagrams, 6 X 9 
in. New York, Academic 
1961. Price, $14.50. 


Press, Inc., 


Twenty-nine papers are included in these 


proceedings of a symposium held in March, 


1960. Although the majority of the contribu- 
tors are from the United States, there are 
papers from England, Germany and Sweden. 
Emphasizing the rapid rate of advancement 
in this field, the papers deal not with the 
question of whether free radicals occur in 
biological systems, but with specific mecha- 
nisms in which they do occur (for example, 
enzyme-substrate compounds, photodynamic 
systems, cytochrome, peptides, etc.). 


Book Nores 


SEQUENTIAL DecopinG, by John M. Wozen- 
craft and Barney Reiffen. 67 pages, dia- 
grams, 5} X 9 in. New York, John Wiley 
& Sons, Inc.; Cambridge, The Technology 
Press; 1961. Price, $3.75. 


The six chapters in this monograph present 
a unified approach to the communications 
problem of coding and decoding, from a proba- 
bilistic point of view. Chapter 1 relates 
coding theory to data 
Chapter 2, on block codes, surveys prior re- 
sults essential to the sequel. Sequential de- 
coding and convolutional encoding are dis- 
Experimental 


communication. 


cussed in Chapters 3 and 4. 
results and modifications of the decoding 
procedure are contained in Chapter 5 on 
“Simulation,” and the final chapter considers 
applications of the concepts to channel models 
and possible engineering applications. This 
is not a text, but a research monograph. 


TRANSMISSION OF INFORMATION: A STATISTI- 
CAL THEORY OF COMMUNICATIONS, by 
Robert M. Fano. 340 pages, diagrams, 
6 X9in. New York, John Wiley & Sons, 
Inc.; Cambridge, The Technology Press; 
1961. Price, $7.50. 


Specifically directed to graduate students 
and engineers interested in electrical com- 
munications, this book is the most up-to-date 
one available in the field of information theory. 
The nine chapters cover transmission of in- 
formation, a measure of information, simple 
message ensembles, discrete stochastic sources, 
transmission channels, channel encoding and 
encoding for binary symmetric 
distributions, and 
constant channels. 
problems for all 


decoding, 
multinomials 
discrete, 


channels, 
encoding for 
Appendix A 
the chapters; 
pages of entropy tables; and 
table of the Gaussian 


contains 
Appendix B_ presents 15 
Appendix C 
is a distribution 
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PUBLICATIONS RECEIVED 


[INDUSTRIAL FURNACES, VOLUME 1, by W. Trinks and Matthew Mawhinney. Fifth edition, 
revised, 486 pages, diagrams, 6 X 9 in. New York, John Wiley & Sons, Inc., 1961. 
Price, $17.00. 


Toor STEEL SIMPLIFIED, by Frank R. Palmer and George V. Luerssen. Third edition, com- 
pletely revised, 595 pages, illustrations, 6 X 9 in. Reading (Penna.), The Carpenter 
Steel Company, 1960. Price, $2.50 (postpaid in U. S.). 


STRUCTURAL ENGINEERING FOR PROFESSIONAL ENGINEERS’ EXAMINATIONS, by Max Kurtz. 
341 pages, diagrams, 54 X 8in. New York, McGraw-Hill Book Co., Inc., 1961. Price, 
$9.00. 


Marinas, by Charles A. Chaney. Second edition, 247 pages, illustrations, 8} X 11 in. 
New York, National Association of Engine and Boat Manufacturers, Inc., 1961. No 
price. 


INTRODUCTION TO ENGINEERING MECHANICS, by John V. Huddleston. 474 pages, diagrams, 
6 X9in. Reading (Mass.), Addison-Wesley Publishing Co., 1961. Price, $9.75. 


THe Impact OF THE NEW Puysics, by Frank Hinman. 168 pages, diagrams, 5} X 8} in. 
New York, Philosophical Library, Inc., 1961. Price, $4.50. 


WATER TREATMENT FOR INDUSTRIAL AND OTHER Uses, by Eskel Nordell. Second edition, 
revised, 598 pages, illustrations, 6 X 9 in. New York, Reinhold Publishing Corp., 1961. 
Price, $12.00. 


PILE FounpatTions, by Robert D. Chellis. Second edition, revised, 704 pages, illustrations, 
diagrams, 6 X 9in. New York, McGraw-Hill Book Co., Inc., 1961. Price, $16.00. 


Second edition, revised, 168 pages, 
Price, $5.50. 


TECHNICAL REPORT WRITING, by Fred H. Rhodes. 
54 X 8} in. New York, McGraw-Hill Book Co., Inc., 1961. 


A First CoursE IN MATHEMATICAL STATISTICS, by C. E. Weatherburn. Second edition 
reprinted with corrections, 277 pages, 54 X 8} in. New York, Cambridge University 
Press, 1961. Price, $2.75 (paper). 


Tue METHODS OF PLANE PROJECTIVE GEOMETRY, by E. A. Maxwell. Reprinted, 230 pages, 
54 X 84in. New York, Cambridge University Press, 1961. Price, $1.95 (paper). 


Rapio ContROL MANUAL, by Edward L. Safford, Jr. 192 pages, illustrations, 54 X 8} in. 
New York, Gernsback Library, Inc., 1961. Price, $3.20 (paper). 
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Television Weather Satellite May 
Spot Hurricanes.—A new chapter in 
one of the nation’s most successful 
space programs opened with the 
launching of the third Tiros ‘‘weather 
eye” satellite to observe weather pat- 
terns and changes around the world 
and especially to report on the origin 
and development of hurricanes in 
the Atlantic and Caribbean areas. 

The newest weather satellite was 
built for NASA’s Goddard Space 
Flight Center by the same Radio 
Corporation of America engineering 
team that provided the Tiros | and 
II satellites and their associated 
ground systems. Where the first two 
models carried one wide-angle tele- 
vision camera for large-area viewing 
and one narrow-angle unit for de- 
tailed cloud study, the new Tiros 
carries two wide-angle cameras. It 
also differs from its predecessors in 
carrying a new experimental infrared 
sensing system developed by Pro- 
fessor Vernon Suomi, of the University 
of Wisconsin, to measure the heat 
radiation coming earthward from the 
sun, and the amounts absorbed and 
reflected by the earth and the at- 
mosphere. 


New Features on Tiros III 


Identical in size and shape to the 
first two Tiros satellites, the 285-lb. 
Tiros III differs in only a few details 
from the earlier types. In addition to 
the substitution of a second wide- 
angle camera for the earlier narrow- 
angle type, and the addition of the 
new infrared experiment, it carries 
improved remote control programmers 
for operating the cameras and other 
electronic equipment in the satellite. 
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Each of the two cameras will photo- 
graph cloud patterns over an area 
about 700 miles on a side for trans- 
mission to the ground stations. Ex- 
cept for the replacement of the previ- 
ous narrow-angle lens by a wide-angle 
type, the duplicate television camera, 
tape recording and transmission sys- 
tems in the new satellite are identical 
to those in the earlier Tiros models. 

The new infrared equipment carried 
for the first time on Tiros III consists 
of an omnidirectional sensor system 
carried on outriggers extending from 
opposite sides of the satellite’s circular 
base. The sensors will determine the 
thermal radiation characteristics of 
the earth by measuring the heat radi- 
ated toward earth by the sun, and 
the amount that is re-radiated from 
the earth’s surface and atmosphere. 
In addition, Tiros II] carried other 
NASA infrared sensing equipment 
identical to that employed for the 
first time in Tiros II to measure heat 
variations on the earth’s surface and 
in the atmosphere. 

A less obvious change in Tiros III 
is the use of new transistorized circuits 
in the two electronic ‘‘clocks’”’ which 
trigger the operation of the television 
cameras and other equipment at any 
desired point in orbit, in accordance 
with the program that is sent to the 
satellite as it passes within range of 
a ground station. The new circuits, 
replacing an arrangement of moving 
parts in the earlier Tiros clocks, are 
expected to add further to the demon- 
strated reliability of the satellite 
system. 

Two significant changes also have 
been made in the RCA-operated 
ground station complex for Tiros ITI. 
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One of the two main ground stations, 
formerly at Fort Monmouth, N. J., 
has been moved to the NASA Wallops 
Station at Wallops Island, Virginia. 
The second main station remains at 
the Pacific Missile Range in southern 
California. Both are backed up, as 
before, by an additional ground sta- 
tion at the RCA Space Center, 
Princeton, N. J. 

The new satellite is scheduled to 
operate for at least three month in an 
orbit similar to those of Tiros I and 
Tiros Il—about 450 miles above the 
earth, in a path extending north to 
the latitude of mid-Canada, and south 
to the latitude of lower Argentina. 
Within this belt, the satellite will pass 
repeatedly over all parts of the earth 
in successive orbits, televising cloud 
formations and measuring infrared 
radiation. 


Quick-Disconnect Coupling for Liq- 
uid Oxygen and Nitrogen.—Cryo- 
genics, Inc., of Stafford, Va., has 
announced the development of a quick- 
disconnect coupling for use with liquid 
oxygen and nitrogen. The device fills 
the need for a coupling which is simple 
to connect and disconnect, and is 
ruggedized for missile and ground 
handling. 

The oxygen and nitrogen couplings 
are keyed differently so that the pos- 
sibility of erroneous connection is elimi- 
nated. The coupling is manufactured 
from standard hex stock so that no 
special tools are required for assembly. 

The coupling consists of a male and 
female tapered connector which, when 
joined, form a leak-proof seal. The 
heart of the coupling is a yoke-type 
lever. Utilizing a leverage action with 
a high mechanical advantage, the 
coupling can be disengaged instantly 
with a minimum of force. The lever 
pivots on the female coupling and is 
fitted with slots that engage pins on 
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the male assembly. The slots are 
engineered such that by depressing 
the lever, the male coupling is brought 
into position and firmly locked. In 
addition, the lever and the female 
connector are engineered to provide 
a safety lock in the form of a ball 
detent arrangement. This safety fea- 
ture acts to prevent accidental move- 
ment of the lever. 

In the locking operation, the male 
coupling firmly seats against two float- 
ing O-Rings within the female as- 
sembly creating a positive leak-tight 
seal. 

By raising the lever, a cam action 
is produced which disengages the pins, 
resulting in a positive disconnect ac- 
tion even in the face of frost buildup. 

The design of this coupling is such 
that icing of outside surfaces is mini- 
mized. The unit is constructed of a 
low thermal conductivity stainless 
steel. The two O-Rings which act as 
seals minimize heat transfer in that 
they permit a vapor buildup be- 
tween the inside and outside surfaces 
and thereby provide a poor thermal 
conductive path. 

The inner diameters of the male and 
female assemblies are identical and 
there are no obstructions in the flow 
stream, thereby, eliminating pressure 
drops. Couplings with inner diam- 
eters as small as } in. and as large as 
8 in. are available. 

The coupling is suitable for remote 
operation and can also be provided 
with integral check valves to prevent 
fluid flow during disconnect operations. 


Low Temperature Transformer.— 
A new kind of transformer that may 
become a major innovation in power 
technology was described by its in- 
ventor, Dr. Richard McFee, of Syra- 
cuse University and Arthur D. Little, 
Inc., at a meeting of the American 
Institute of Electrical Engineers in 
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Ithaca, New York. The new trans- 
former was developed by Dr. McFee 
at A. D. Little, where it has been 
successfully operated at a level of 
15 kilowatts. 

According to McFee, transformers 
using the new design could be con- 
siderably lighter and less costly than 
conventional transformers. And they 
would not be subject to the substantial 
power losses caused by electrical heat- 
ing, or to the temperature limitations 
of present-day electrical insulation. 

The new design makes use of super- 
conductivity—the ability of certain 
metals to conduct an electric current 
without any resistance at tempera- 
tures near absolute zero (—459.6°F.). 
Superconducting transformers have 
been considered before, but the fact 
that a moderately strong magnetic 
field can quench the superconducting 
state by restoring electrical resistance 
presented an apparently insurmount- 
able obstacle. When an electric cur- 
rent flows through a coil, such as the 
windings of a transformer, powerful 
magnetic fields are generated about 
the windings and quench the super- 
conductivity when they reach a critical 
field strength. 

Dr. McFee discovered that inter- 
leaving layers of the primary and 
secondary windings would keep the 
total magnetic field strength below 
the critical level. Since the current 
in adjacent layers flows in opposite 
directions, the magnetic fields gener- 
ated about these layers are also op- 
posing. In effect, the individual mag- 
netic fields nearly cancel one another 
and eliminate the problem of magnetic 
quenching. 

The Gifford-McMahon helium re- 
frigerator invented at A. D. Little 
provided a solution to another major 
problem in the development of a 
practical superconducting transformer. 
There had not been a simple, reliable 
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method of maintaining the supercon- 
ducting coils at the necessary low 
temperatures for an extended period. 
Now in its final stages of development, 
the new refrigerator will produce tem- 
peratures to —452°F., and will main- 
tain them unattended for months of 
continuous operation. 

Since the cost of refrigerating the 
whole transformer would be excessive, 
only the superconducting coils them- 
selves are cooled, while the cores are 
kept at room temperature. Heat 
leaks into the system are controlled 
by vacuum insulation, liquid-nitrogen 
heat shields, and an ingenious arrange- 
ment of the superconducting input leads. 

Superconducting transformers could 
eliminate a number of major problems 
that face designers of conventional 
transformers. The lack of electrical 
resistance in superconducting metals 
means that little heat is produced. 
Expensive devices for dissipating the 
heat generated in conventional trans- 
formers are therefore unnecessary in 
the new development. For the same 
reason, present power losses from 
heating in the windings would no 
longer exist, and engineers would not 
have to consider the deteriorating 
effect of heat on insulation. 

The problem of excessive voltage 
loads can be handled in the same 
manner as in conventional trans- 
formers. Excessive current loads, 
however, require more elaborate safe- 
guards in the new device because they 
could quench superconductivity and 
put the transformer out of service 
for a minute or more. Such outages 
can be prevented through the use of 
very fast-acting circuit breakers or 
external inductance coils. 


Tiny Experimental Thin-Film Tran- 
sistor.—An ultraminiature experimen- 
tal transistor, so small that as many as 
20,000 can fit on a postage stamp, 
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has been announced by the Radio 
Corporation of America. The tran- 
sistor is made by depositing thin films 
by evaporation on an insulating base 
and is capable of shrinking the basic 
circuitry of a computer to the size 
of a book page. ‘The basic circuitry 
of present-day computers ranges in 
size from the equivalent of a large 
hat box to a walk-in clothes closet. 

Dr. James Hillier, Vice President, 
RCA Laboratories, described the de- 
velopment as ‘‘an important break- 
through in basic electronic circuitry 
that may well have a revolutionary 
impact upon the electronics industry.” 

According to RCA scientists, the 
new device, which has been tested 
successfully at RCA _ Laboratories, 
may open the way to new ultraminia- 
ture mass-production transistor cir- 
cuits for many applications, especially 
in electronic computers and perhaps 
ultimately in other equipment such 
as thin-screen wall-type television 
receivers. 

The active material used in the 
transistor is cadmium sulfide, a com- 
pound with considerably greater in- 
sulating properties than the germa- 
nium, silicon, and other semiconductor 
materials used in standard transistors. 

In making the thin-film transistors, 
an evaporation process is used to de- 
posit successive thin layers of cad- 
mium sulfide and metal on a glass 
plate, creating a device that is only a 
few ten-thousandths of an inch thick. 
In the evaporation process, the cad- 
mium sulfide crystals and the metal 
are heated in successive steps in a 
vacuum, turning to vapor that is 
collected by condensation on the glass 
plate in the same manner as steam 
condenses in a film on a cooler surface 
held over boiling water. 

By using a special mask to cover 
portions of the plate during the proc- 
ess, the metal layers are deposited in 
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a pattern that forms the electrical 
contacts needed to operate the tran- 
sistor. The masking process also can 
be used to produce various patterns 
of connections among many transis- 
tors to complete a desired circuit at 
the same time that the transistors 
themselves are being made. 

The completed transistor is not only 
very tiny, but it also incorporates an 
important operating feature not now 
used in commercial transistors. In 
conventional transistors having com- 
parable functions, electrons flow more 
or less freely through the semicon- 
ductor material between two of the 
contacts, and the third element pro- 
vides control by reducing the flow in 
varying degrees. The operating prin- 
ciple of the experimental thin-film 
transistor is exactly opposite. The 
insulating properties of the cadmium 
sulfide hamper the flow of electrons 
between two electrodes, and the third 
element provides control by increasing 
the flow in varying degrees. 


Condensation Nuclei Detector.— 
General Electric Company has devel- 
oped a device which promises impor- 
tant help in the battle against air 
pollution. It is so sensitive for detect- 
ing tiny airborne particles that it could 
find a single speck lost among 1000 
trillion other specks. This would be 
like spotting a solitary grain of white 
sand on a beach full of yellow grains. 

This “condensation nuclei detector” 
also has important industrial applica- 
tions, although it is not yet com- 
mercially available. ‘Condensation 
nuclei’ are the smallest airborne par- 
ticles in existence. Invisible but ex- 
tremely numerous, they are consid- 
ered a major factor in smog formation. 

The detector could be used to help 
forecast weather conditions by spot- 
ting nuclei present in our atmosphere 
from places hundreds or even thou- 
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sands of miles away. The nuclei 
would move as the air mass moved. 
With the use of auxiliary equipment, 
radioactive particles could be detected 
and, in most cases, traced to their 
source. Monitoring wilderness areas 
for forest fires is another possible use. 
Another unusual application is pre- 
dicting the useful life of certain elec- 
trical equipment on the basis of a 
single operation. An example is that 
of the contact points in an electrical 
relay. Their life could be accurately 
predicted with the points coming 
together just once. There would be 
no visible signs of wear, even under a 
microscope. 

This technique promises an enor- 
mous time saving, for it can accom- 
plish in minutes or hours’ what 
normally takes months. 

The detector works by causing drop- 
lets to be formed artificially around 
each tiny nucleus. The nuclei then 
grow to microscopic size. ‘This is 
large enough to cause scattering of a 
light beam directed into their midst. 

Measured electronically, the light 
scattering indicates the concentration 
of nuclei in the man-made fog. It 
takes but 5-thousandths of a second 
to form a measurable cloud. 

Ability to detect the nuclei with a 
sensitivity of better than one part per 
1000 trillion will prove a boon to all 
forms of air pollution work and should 
make possible the design of air filters 
with an effectiveness far exceeding any 
now in existence. 

Monitoring air pollution, including 
the exhaust from automobiles, would 
be another important job for a nuclei 
detector of this high accuracy. It 
has been estimated that air pollution 
is presently costing the American 
public more than $10 billion a year in 
corrosion, crop damage and impair- 
ment of health. 

Continuous operation of the device 
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makes it of particular value in detect- 
ing and identifying gases. The gases 
cannot be detected directly, but an 
indirect measurement is possible if the 
gas molecules are combined with sub- 
stances to give them a solid or liquid 
form. This is accomplished through 
the use of “‘converters’’ especially de- 
signed for the requirements of each 
basic group of gases. 

The extreme sensitivity of the de- 
tector is due in large measure to its 
continuous operation, which also per- 
mits its use for automatic monitoring 
of the air in a remote location. The 
continuous operation is achieved by 
forming and measuring five to eight 
cloud masses a second inside the in- 
strument’s ‘“‘cloud chamber.”’ Previ- 
ous laboratory models took inter- 
mittent measurements of the air, 
spaced more than a minute apart, and 
yielding results that had to be averaged 
out. 


Micrometeorite Layer Discovered. 
—The Air Force has announced the 
discovery of what appears to be a 
dense band of micrometeorites which 
envelopes the earth at extreme 
altitudes. 

This discovery was made recently 
by the Air Force Cambridge Research 
Laboratories, Laurence G. Hanscom 
Field, Bedford, Mass., when the con- 
tents of a novel rocket nosecone 
specially designed for trapping micro- 
meteorites were examined. The nose- 
cone was launched June 6 at White 
Sands, New Mexico. The density 
of the micrometeorites in this layer 
was totally unexpected. About 10 
of these small particles struck each 
square centimeter of the detecting 
surfaces each second. 

The dust particles, measured in 
microns (thousandths of a millimeter) 
apparently exist as a band about the 
earth that is independent of latitude 
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or longitude, Dr. R. K. Soberman, 
AFCRL project scientist, said. The 
precise altitude at which these micro- 
meteorites are concentrated will be 
determined by future tests now planned 
by AFCRL. 

The special nosecone for trapping 
and recovering these micrometeorites 
was designed by scientists and engi- 
neers of AFCRL and the Aerojet 
General Corp. It consisted of a pod- 
like arrangement made up of indi- 
vidual petals or leaves. At an altitude 
of 47 miles these leaves opened like 
a blossoming flower. 

The nosecone was carried to a maxi- 
mum altitude of 102 miles over the 
desert at White Sands, New Mexico. 
The leaves closed again when the 
rocket descended to 65 miles. Alto- 
gether the nosecone remained opened 
for about 4 minutes. 

Two types of micrometeorite detec- 
tors were used. ‘The first type con- 
sisted of a triple layer of physically 
separated Mylar and Plexiglas sheets. 
Each sheet measured about 2” X 3”. 
Detectors of this type were located 
on each of the leaves of the nosecone. 
The top layer consisted of a Mylar- 
film } mil (1/4000ths in.) thick. The 
second layer, also of Mylar was 1 mil. 
The bottom layer was a }” thick sheet 
of Plexiglas. The micrometeorites 
traveling at speeds up to 47 miles a 
second, passed through the 2 Mylar- 
film layers, creating many holes. 
Most of these were microscopic in size, 
although many were visible to the 
naked eye. When they struck the 
harder Plexiglas layer below, small 
craters, also readily visible, were 
created. 

In addition to detectors of this type, 
relatively thick films of 3 plastics of 
varying hardness—Millipore, Form- 
var, Lucite—were used. Craters simi- 
lar to those formed in the Plexiglas 
were detected in the Lucite film. 
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Most of the evidence of micromete- 
orite contact, however, was apparent 
only through microscopic examination. 

Few micrometeorites themselves 
were collected since, for the most part, 
they vaporized on contact with the 
detecting surfaces. Some of these 
extra-terrestrial particles, however, 
were collected intact. Of particular 
interest to AFCRL scientists is the 
residual bits of meteorite material 
lining the walls of the craters. These 
residual materials will be analyzed 
by physical, chemical and radioactive 
techniques. 

Dr. Soberman suggested that the 
micrometeorite layer is probably 
formed by electrostatic trapping. 
Complete analysis, Dr. Soberman said, 
may take several years. The analysis 
will be done for the most part by 
electron microscopy techniques. The 
electron microscope must be used 
because many of the individual par- 
ticles to be examined are smaller than 
the wave length of light. Complete 
examination of the 1 square inch of 
material, for example, may take a 
single scientist one full year. A large 
number of samples are available and 
AFCRL will furnish these to inter- 
ested and qualified scientists. 

One limitation on the present study 
is that there is no way of knowing 
with precision the altitude at which 
the meteorite layer is concentrated. 

AFCRL hopes to determine this 
when a specially instrumented rocket 
is launched from Eglin AFB, Fla., in 
the near future. Improved micro- 
phone detection techniques which 
record the impact of these tiny par- 
ticles on special metallic surfaces will 
be used in this test. 


Instrumentation Rack for Sounding 
Rockets.—The Geophysics Corpora- 
tion of America, Bedford, Mass., has 
developed a standardized instrumen- 
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tation rack for upper atmosphere 
sounding rockets. The rack 
signed to accommodate equipment 
payloads of variable size and configura- 
tion and is able to withstand severe 
conditions of shock and vibration. 

Built to package instrumentation 
payloads for rockets launched both 
in the United States and in foreign 
countries, the new rack is made of 
6061-T6 aluminum with stainless steel 
hardware. It has four struts, one of 
which serves as a wiring channel and 
contains a series of spaced openings 
for plug-in assemblies. Although its 
length may be varied according to 
the needs of each experiment, the 
model available for ASP, CAJUN or 
APACHE rockets stands 25,'s inches 
high and weighs 8.5 lb. empty. 

Evidence of the new rack’s sturdi- 
ness was demonstrated recently in a 
rocket launching in which the rocket 
blew up 25.6 sec. after launch. Com- 
plete telemetry recordings continued 
to be received from the instru- 
mentation aboard until the equipment 
impacted approximately 297.4 sec. 
later. From the telemetry records 
it was calculated that the rack and 
equipment was subjected to 1,000 
“G”’ forces of radial acceleration, and 
to an unmeasured but extremely high 
amount of vibration 


is de- 


Multiple Diode.—The Radio Cor- 
poration of America recently intro- 
duced the electronic industry's first 
multiple germanium switching diode 
in a new bid to simplify the manu- 
facture of computers. These new 
multiple diodes consist of two and 
three semiconductor diodes of ad- 


vanced technology, packaged as single, 
integrated units, according to Dr. 
Alan M. Glover, Vice President and 
General Manager, RCA Semiconduc- 
tor & Materials Division. 

He said they are expected to sim- 
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plify the assembly of computers by 
reducing the number of separate com- 
ponents and electronic connections 
required to build them. 

Dr. Glover explained that diodes 
are small, individual, two-terminal 
devices used to help translate into 
computer electronics the principle that 
every question, no matter how com- 
plex, can be answered by reducing 
it to a series of “‘yeses’’ and “noes.” 

In conjunction with transistors, he 
said, diodes perform these “‘yes’”’ and 
“no” functions by opening or closing 
thousands of “logic gates” (circuits) 
that control the flow of coded elec- 
tronic signals through the computer. 

“Building this logic network can be 
a very complicated time-con- 
suming job,”’ Dr. Glover said. ‘‘How- 
ever, if individual diodes are replaced 
by multiple diode assemblies, the job 
can be greatly simplified while giving 
rise to such other potential benefits 
as quicker assembly schedules, greater 
space savings and lower production 
costs.” 

Commenting on the advanced tech- 
nology embodied in the new diodes, 
Dr. Glover disclosed that they are 
manufactured by a _ new _ process, 
unique to the diode field, that provides 
very high switching speeds as well 
as great mechanical and electrical 
stability. 

“The simplicity and reliability of 
this unique process have been matched 
with a simplicity and reliability of 
design which,” he said, “gives our 
new germanium diodes wide utility 
with high performance and low cost. 
They should prove of exceptional 
value to the entire computer industry.” 


Automated Production of Deposited 
Carbon Resistors.——A wholly new 
concept in the automatic manufacture 
of electronic components that com- 
bines statistical quality control with 
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fully mechanized processing—all un- 
der the guidance of a digital computer 
—was given its first public demon- 
stration recently at the North Caro- 
lina Works of the Western Electric 
Company. 

The block-long array of intercon- 
nected robots, observed in operation 
by newspaper reporters and technical 
editors, manufactures deposited car- 
bon resistors. The resistors are elec- 
tronic components that must meet 
severe electrical requirements, yet 
operate thousands of years on the 
average without failure. So great is 
the reliability demanded of the prod- 
uct that manual techniques of manu- 
facture are impractical. ‘The resistors 
are used in immense quantities and in 
a variety of electronic equipment in- 
cluding radar and missile control 
systems. 

In addition to controlling step-by- 
step details of manufacture, the com- 
puter can accept a production schedule 
for a full month’s output and there- 
after automatically issue detailed in- 
structions to the fabricating machines 
to make and package any or all of 
four basic power sizes of resistors in 
such quantities and resistance values 
as may be desired. 

According to engineers who devel- 
oped the automated line, it is ‘the 
first time that a computer has been 
used in conjunction with a statis- 
tical quality control system to con- 
trol a series of fabricating machines 
which produce precise, electronic com- 
ponents.” 

The resistors look like tiny black 
firecrackers with copper wire “‘fuses”’ 
extending from the ends. Each con- 
sists of a short ceramic rod or core 
which is first coated with carbon. A 
conductive gold coating is applied to 
each end of the core for attaching a 
cap and wire terminal. Then a 
spiralled groove is cut into the carbon 
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film, barber-pole fashion, to change 
the electrical path on the core and 
raise the resistance to the desired 
level. Finally, the resistor is encased 
in a cylindrical, epoxy plastic cover 
for protection. The machines work 
on a three second cycle, producing 
a resistor during each cycle. Maxi- 
mum production rate is 1200 finished 
resistors per hour. 

There are eleven machine stations 
on the line, in addition to the com- 
puter and control cabinets. Three of 
these stations are used only for testing. 
Length of the line is 110 feet, with a 
maximum machine width of about six 
feet. 

To the observer's eye, the line’s 
eleven machine stations have the trim, 
console look characteristic of modern 
industrial machines. One, in par- 
ticular, puts on a brilliant show. 
Known as the “‘sputtering machine,” 
it sprays the ends of the resistors with 
microscopic particles of pure gold 
inside gas-filled bell jars. The dis- 
charge of high-voltage electricity that 
propels the precious metal creates a 
veritable aurora borealis in miniature 
of flickering reds and purples. 

In missile systems, deposited carbon 
resistors are required by the millions 
and they must work under extreme 
conditions of cold, heat, humidity, 
vibration and shock. In certain de- 
fense equipment, a failure rate of no 
more than one per 200 million hours 
(some 23,000 years) of operation is 
permissible. 

It is mainly this reliability require- 
ment which triggered the design and 
development of the machines shown 
here today. Until now, individual, 


precision products have been made 
by manual or semi-automatic proc- 
esses that invited contamination from 
handling and related shortcomings of 
human control. 
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By adopting advanced automation, cannot match, advantages that are 
contamination disappears as a prob- reflected in low unit costs. 
lem, production is vastly increased Western Electric will continue the 
and, most vital, a level of reliability use of manual methods for the manu- 
is attained that manual methods facture of conventional resistors. 


YOU CAN ADVANCE SCIENCE EDUCATION 


Today, more than ever before in its 138-year history, there is vital need for 
The Franklin Institute effectively to promote education in science and technology. 
It is imperative that we meet this challenge by providing adequate educational op- 
portunities in these fields. This requires vision, objective planning, and money. 
We have more than enough of the first two requisites, but far too little of the third. 

Our programs are aimed at professional scientists and industry, as well as the lay 
public and young people seeking inspiration and guidance in choosing a career. 
The Institute’s educational programs are impressive. They begin with students in 
the early grades of our elementary schools and continue throughout an individual's 
professional or industrial life. With more funds at our disposal, the scope and 
vigor of these activities could be greatly increased and increasingly effective. 

The Franklin Institute is not richly endowed. It is a non-profit organization, 
depending for encouragement and support on an understanding public. Capable 
and conservative management assures wise administration of all funds. 

Your gift or bequest, large or small, will be deeply appreciated and will be used 
effectively to broaden the Institute’s educational usefulness. There is a warm 
satisfaction in giving financial support to an organization that has pioneered in, and 
is dedicated to, the advancement of science and technology. 

When property is transferred, title should be in the name of The Franklin 
Institute of the State of Pennsylvania for the Promotion of the Mechanic Arts. 

The Secretary of The Franklin Institute will gladly furnish you with additional 
information. Write to him at The Franklin Institute, Benjamin Franklin Parkway 
at Twentieth Street, Philadelphia 3, Pennsylvania. 

If you wish to give securities: (1) Send the stock certificates to The Franklin 
Institute, Philadelphia 3, Pennsylvania, and leave the assignment space blank. (2) 
Execute a stock power and mail this to the Institute separately, leaving the assign- 
ment space blank. (The Institute, upon request, will gladly furnish blank stock 
power forms.) 
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STAFF POSITIONS 


. are available for qualified scientists and engi- 
neers in the Laboratories of The Franklin Insti- 
tute. These Laboratories comprise one of the 
country’s finest independent research organiza- 
tions, and provide R&D services to both in- 
dustry and government, in a variety of fields: 


chemistry, solid state sciences, applied physics, 
operations analysis, engineering psychology, 
and chemical, electrical, mechanical and nu- 
clear engineering. 


There are openings, in groups working on chal- 
lenging problems in both research and develop- 
ment, for men with proven abilities in any of 
these fields. 


THE FRANKLIN INSTITUTE LABORATORIES 
FOR RESEARCH AND DEVELOPMENT 


20th Street and Benjamin Franklin Parkway 
Philadelphia 3, Pennsylvania . .. LOcust 4-3600 


Send complete resumé to 
Mr. John E. Christ, Director of Personnel 


Or write for a brochure describing activities 
of the Laboratories 
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BUSINESS 
MATH 
ENGINEERING 


Now in the fourth year of operation, The Franklin Institute Com- 
puting Center has shared its technical know-how, modern electronic 
equipment and experienced personnel with hundreds of progressive 
industries and government agencies across the nation. 


We offer the services of creative people, highly skilled in their respec- 
tive fields and ably trained in the application of these skills to the ever 
expanding area of electronic computers and data processing systems. 


This staff is now available for analysis, system design, programming 
or coding of projects of unlimited scope or context. Input to our large 
scale computer and peripheral equipment is acceptable in any form. 
Results are provided on cards, plastic or metallic tape, and in com- 
pletely edited printed copy. Our extensive library of automatic coding, 
engineering, data processing and mathematical routines is available to 
all users, and machine time is provided with or without the services 
of programming personnel. 


THE FRANKLIN INSTITUTE 
Computing Center 


20th Street and Benjamin Franklin Parkway 
Philadelphia 3, Pa. LOcust 4-3600 
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